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1. Introduction 

The characterization of states of matter in terms of 
spontaneously broken symmetries and the corresponding 
order parameters, known as Ginzburg-Landau theory, is 
one of the most important concepts in condensed matter 
theory. Perhaps the best known and studied example 
is that of magnetic order, for instance originating from 
strong electronic interactions. Remarkably, there also 
exist insulating phases of matter which do not break 
any symmetries, yet are distinct from simple band 
insulators. The first such state to be discovered was 
the integer quantum Hall state of a two-dimensional 
electron gas in a strong magnetic field [1], see [2] for 
an introduction. Its quantized Hall conductivity (In two 
dimensions, conductance and conductivity are identical.) 
can be explained by the existence of a topological invariant 
(the Chern number) calculated from a single-particle 
Hamiltonian [3-5]. Because this number cannot change 
under adiabatic transformations of the Hamiltonian, the 
integer quantum Hall state is topologically different from 
a trivial band insulator. In particular, a transition between 
these two states can only occur via a closing of the band 
gap. As a result, metallic states emerge at the edge 
of quantum Hall samples [6]. Moreover, while atomic 
insulators may be described as a product state of completely 
local electronic wave functions, topological states show 
a distinct response to boundary conditions [7]. Haldane 
[8] showed that a quantum Hall state with broken time- 
reversal symmetry can not only be achieved with an external 
magnetic field, but also with a periodic magnetic flux. 

Topological insulators represent another class of 
states that are topologically distinct from simple band 
insulators. They are symmetric under time reversal, and 
can therefore be realized experimentally without magnetic 
fields. Topological insulators were first discovered in 
a theoretical model for graphene [9, 10], in the context 
of the search for an intrinsic quantum spin Hall effect. 
Kane and Mele [9, 10] showed that the intrinsic spin- 
orbit coupling gives rise to an insulating state with 
a topological band structure, and a nonzero spin Hall 
response related to metallic edge states that are protected 
against perturbations as long as time-reversal symmetry is 
not broken. Topological insulators in three dimensions were 
predicted soon after [11], and are by now well confirmed by 
experiments [12]. Two-dimensional topological insulators, 
or quantum spin Hall insulators, have close conceptual 
relations to integer quantum Hall states: in the simplest 
case, the theoretical model of Kane and Mele [9, 10] 
is completely equivalent to two decoupled copies (with 
opposite chirality for the two spin directions) of the Haldane 
model for the integer quantum Hall effect [9]. 

While the spin-orbit gap in graphene is too small 
to permit experimental observation (about 10 -3 meV, see 
for example [13]), the quantum spin Hall state has been 
observed in CdTe/HgTe/CdTe heterostructures [14-16]. As 



first predicted theoretically [17], the low-energy physics of 
this system is determined by the gap at the T point. If 
the HgTe quantum well has a thickness below a critical 
value d c , the band structure is dominated by the surrounding 
CdTe layers, and the normal regime, with the s (Tq) bands 
lying energetically above the p (r 8 ) bands, is realized. If 
the thickness exceeds the critical value, the band structure 
becomes dominated by that of the quantum well, and band 
inversion takes place. Bernevig et al [17] predicted samples 
with such an inverted band structure to be quantum spin 
Hall insulators with a band gap of tens of meV, which was 
soon verified experimentally [14, 15]. 

In the absence of interactions, topological insulators 
can be understood in terms of single-particle Hamiltonians, 
allowing for a full classification based on time-reversal, 
particle-hole, and chiral symmetry [18-20]. Recently, an 
extension to include crystal symmetries has been given 
[21]. This review is primarily concerned with correlation 
effects, which have been explored in the context of 
topological insulators in several different settings. The 
most direct way to theoretically explore correlation effects 
is to consider the impact of electron-electron interactions 
on a noninteracting quantum spin Hall state by adding 
interactions to models with intrinsic spin-orbit coupling. 
Remarkably, some of these models permit the application 
of exact numerical methods. With regard to experiment, 
such studies are relevant for materials with significant spin- 
orbit and electron-electron interactions. This situation is 
expected to occur in, for example, the transition-metal 
oxide Na 2 Ir0 3 [22]. The interplay of strong correlations 
and spin-orbit coupling in this and other materials may 
allow for realizations of complex spin models [23-26] and 
hence may provide a route to quantum spin liquids [27]. 

Even in materials with rather weak intrinsic spin-orbit 
coupling, robust quantum spin Hall states can emerge from 
strong electronic correlations via dynamically generated 
spin-orbit coupling as a result of spontaneous symmetry 
breaking [28]. The resulting topological Mott insulators, 
states with a band gap generated by interactions and edge 
states protected by time-reversal symmetry, have so far been 
explored only at the mean-field level. The term topological 
Mott insulator is also used in the literature to refer to 
exotic states in which the interplay of spin-orbit coupling 
and electronic correlations leads to spin-charge separation, 
gapped charge excitations, and a topological band structure 
carried by low-energy spinon excitations [29]. 

A third direction of research is concerned with 
fractional topological insulators. The fractional quantum 
Hall state is a famous example of a topological phase 
emerging from strong electronic interactions; for a concise 
introduction see [2]. In the latter, strong repulsion 
among electrons occupying partially filled, flat Landau 
levels leads to intriguing insulating states with quasi- 
particle excitations that carry fractional charge and have 
fractional statistics. Fractional quantum Hall states have 
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been a major research topic for several decades, with 
enormous progress in theory and experiment, and potential 
applications in quantum computing [2]. Remarkably, much 
of the physics of the fractional quantum Hall effect can 
in theory also exist in Chern insulators [30], described 
by models that break time-reversal symmetry but have 
zero net magnetic flux, and are hence similar in spirit 
to Haldane's model for the integer quantum Hall effect 
[8]. Proposals for the experimental realization of fractional 
Chern insulators involve materials with flat bands that allow 
for an enhancement of electron-electron interactions [31]. 
In close analogy to the noninteracting case, fractional 
quantum spin Hall insulators can be constructed by, in 
the simplest case, combining two fractional quantum Hall 
states [32]. The theory of time-reversal invariant insulators 
with protected edge states is one of the most active current 
research areas (see section 4.5). 

In the context of correlated electrons, the concept 
of topological order has yet another, important meaning. 
Quite generally, a topologically ordered state is robust with 
respect to adiabatic transformations of the Hamiltonian. 
In the case of the Z 2 topological insulators, as described 
for example by the Kane-Mele model [9, 10], this 
robustness is mathematically expressed by a topological 
invariant. However, even in the presence of electron- 
electron interactions, this state — having its origin in a 
topological band structure — is adiabatically connected to 
a noninteracting band insulator. A classification of such 
states based on Chern-Simons theory has recently been 
given [33]. Topological band insulators are distinct from 
the class of topological states that emerge from interactions 
and are not adiabatically connected to a noninteracting 
state. A particularly fascinating example are quantum spin 
liquids, which are insulating states — arising from quantum 
fluctuations — that do not break any symmetries but have 
highly nontrivial correlations. For the latter, the term 
topological order implies [34] (i) the existence of a gap 
to excited states, (ii) a topological ground-state degeneracy 
in a periodic system (a torus geometry in two dimensions), 
and (iii) the existence of fractional quasi-particle excitations 
with Abelian or non-Abelian statistics. Whereas (i) is also 
typical of topological band insulators, (ii) and (iii) are not. 
On the other hand, fractional Hall liquids and hence also 
fractional topological insulators possess topological order. 
A review of these concepts and the connections between 
topological insulators and quantum spin liquids has been 
given by Fiete et al [35]. 

Whereas a rich variety of materials has emerged 
that support three-dimensional topological phases [12,36, 
37], experiments for the two-dimensional case are so far 
restricted to HgTe quantum wells [14-16], and do not reveal 
any substantial electronic correlation effects. Consequently, 
the study of correlated quantum spin Hall insulators has 
so far been a predominantly theoretical effort. The search 
for candidate condensed matter settings is guided by the 



fact that the spin-orbit coupling increases with the atomic 
number. For example, strong and comparable electron- 
electron and spin-orbit interactions are expected to exist 
in Ir-based transition metal oxides [22], with numerical 
evidence for a strongly correlated quantum spin Hall state 
[22], and the possibility of a topological Mott insulator 
phase [29]. Another useful experimental setting are 
heterostructures of transition-metal oxides, for example a 
bilayer of LaNi0 3 sandwiched between LaA10 3 layers, 
in which spin-orbit coupling is mimicked by interaction- 
induced spontaneous ordering of complex orbitals [38, 
39]. Interaction induced topological phases may also exist 
in stacked graphene [40,41], whereas the intrinsic spin- 
orbit coupling may be sufficiently enhanced in decorated 
graphene [42] or in molecular graphene [43] to make a 
topological phase observable. Although concrete materials 
for fractional topological insulators have not yet been 
identified, several candidate systems for Chern insulators 
are known, see [31] and references therein. Cold atoms 
in optical lattices, or on atom chips, provide an alternative 
route toward correlated topological states. Using suitable 
artificial gauge fields [44], Z 2 topological insulators with 
tunable interactions may be realized [45,46]. Finally, 
correlated topological states may also be engineered with 
photons in cavity arrays, see [48, 49] and references therein. 

Despite the rather short history of topological insula- 
tors, the literature is remarkably rich. Given the background 
of the authors, this review aims at providing an overview of 
work on electronic correlation effects in two dimensions, 
with a focus on models with intrinsic spin-orbit coupling 
and numerical studies. It thereby complements earlier re- 
views by Hasan and Kane [36], Moore [50], and Qi and 
Zhang [12], but omits other fascinating topics such as three- 
dimensional topological insulators [12, 36, 37], topological 
Kondo insulators [51], or topological superconductors [12, 
18,19,36]. For a self-contained presentation, it also in- 
cludes an introduction to noninteracting quantum spin Hall 
insulators using the example of the Kane-Mele model. The 
review is organized as follows. In section 2, the most widely 
used theoretical models are introduced. A discussion of 
noninteracting quantum spin Hall insulators is provided in 
section 3. Correlation effects on bulk properties are con- 
sidered in section 4, whereas correlated helical edge states 
are the topic of section 5. The calculation of topological in- 
variants for correlated systems is reviewed in section 6, and 
section 7 gives conclusions and an outlook. 

2. Theoretical models 

Since the proposal by Kane and Mele [9, 10], a number of 
theoretical models for the quantum Hall and spin Hall effect 
have been discovered. The general structure of such models 
and possible realizations in semiconductors were discussed, 
for example, in [52]. In many cases, spinless quantum Hall 
models can readily be generalized to describe time-reversal 
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invariant quantum spin Hall states by introducing a spin- 
orbit gap with a spin-dependent sign. Here the emphasis 
is on some of the most popular models that support a 
topological state originating from spin-orbit coupling in the 
absence of electronic interactions. Overviews of theoretical 
models can also be found in [7, 35]. 

A common feature of the quantum spin Hall models 
considered here is that the Hamiltonians may be expressed 
in terms of 4 x 4 Dirac matrices [7,9]. Given 
inversion symmetry, such a representation is particularly 
useful to determine the topological invariant in the 
absence of interactions [7,9]. The corresponding matrix 
representations of the KM and the Bernevig-Hughes-Zhang 
(BHZ) model can be found in [7]. To study local, Hubbard- 
type interactions, and also for numerical methods, the tight- 
binding real-space representations are more appropriate. 

Because of the topological character of the quantum 
spin Hall state, the details of the Hamiltonian play a minor 
role, and different models share similar low-energy physics. 
In particular, universal field theory descriptions can be 
obtained using Chern-Simons theory for the gapped bulk 
[53], and helical Luttinger liquid theory for the gapless 
edges [54,55]. The above-mentioned requirement of at 
least four orbitals per unit cell (here the term orbital may 
also refer to spin or pseudospin) or, equivalently, four states 
at a given momentum k in the noninteracting Bloch picture, 
can be understood as follows: the two spin directions are 
related by time-reversal symmetry and hence degenerate. 
A band insulator, with a gap between filled valence bands 
and empty conduction bands and hence not adiabatically 
connected to the vacuum, then requires two sets of Kramers 
degenerate bands, and hence four bands in total [56]. 

While some features are generic for quantum spin Hall 
models, there are also model-specific aspects related to, 
for example, the crystal lattice. An important example 
is the potential quantum spin liquid phase of the KMH 
model, which was first reported for the Hubbard model on 
the honeycomb lattice [57]. This lattice has the smallest 
possible coordination number z = 3 in two dimensions, 
rendering quantum fluctuations particularly important. The 
honeycomb lattice is also at the centre of attention in 
condensed matter physics due to tremendous progress in 
preparing single layers or bilayers of graphene [58, 59]. 

Much of the work reviewed here uses the paradigmatic 
framework of a local Hubbard interaction to study 
electronic correlation effects [60]. This choice, usually 
motivated by simplicity, is assumed to capture the dominant 
physics. Models with nonlocal interactions play a key role 
for topological Mott insulators with dynamically generated 
spin-orbit interactions (see section 4.4), and for spinless 
models of Chern insulators (see section 4.5). Nonlocal 
interactions may also give rise to interesting new physics 
in helical liquids [55]. 




Figure 1. (a) The honeycomb lattice with lattice constant a consists 
of two sublattices A, B and is spanned by the basis vectors a\ = 
^a(3, \/3)> o-2 = 51(3, — %/3) (in the conventions of [59]). Nearest- 
neighbour lattice sites are connected by the vectors <5i = 2<i(l,\/3), 
52 = § a (li — Vo), and £3 = a(~ 1,0). The hexagonal first Brillouin 
zone contains the two nonequivalent Dirac points K 
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2.1. Kane-Mele-Hubbard model 

The KM Hamiltonian was originally derived as a model for 
graphene, and can be written in the form [9, 10] 



H, 



KM 



-t 



(1) 



+ iA 



aZ Jgt (sxd l] ) z c J 

(id) 



Here cj = (c^, c^) is a spinor, and cj CT creates an electron 
in a Wannier state at site i with spin a. Pairs of 
nearest-neighbour and next-nearest-neighbour lattice sites 
are indicated by the symbols and which also 

implicitly include the Hermitian conjugate terms. 

The KM model (1) is defined on the honeycomb 
lattice spanned by the lattice vectors o,, a 2 , as shown 
in figure 1(a) (see caption for details). Nearest-neighbour 
sites are connected by the vectors ±5,, i = 1,2,3. 
The honeycomb lattice has two lattice sites per unit cell, 
and together with the spin degree of freedom, there are 
four orbitals per unit cell. The Brillouin zone, shown in 
figure 1(b), contains the Dirac points K and K'. 

The first term in (1) is the usual tight-binding hopping 
term between sites on different sublattices, which leads 
to the celebrated graphene band dispersion (with lattice 
constant a = 1) [61] 



<k)=±\g k \ 



(2) 



±t 



3 + 2cos(v / 3fcJ 



4 cos (3fc x /2) cos(V3\/2) 



1/2 



with Dirac cones at K and K'. 

The spin-orbit interaction (~ A go ) is described by a 
complex-valued next-nearest-neighbour hopping term with 
a sign ±1 depending on the sublattice, the direction of the 
hop (i.e., left turn or right turn), and the spin orientation, 
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+i\ SO (T Z -i\ SO cr z 

Figure 2. Sign structure of the spin-orbit term in the Kane-Mele model (1). 



existence of a Mott metal-insulator transition at a finite 
critical U on the honeycomb lattice is a result of the 
vanishing density of states at the Fermi level, N(ui) ~ 
\oj\ [59], which renders the logarithmic instability toward 
magnetic order ineffective. This stability with respect 
to perturbations is a general property of Dirac points. 
In contrast, symmetry breaking and topological phases 
can arise for arbitrarily small interactions in models with 
quadratic band crossing points [64]. 



as illustrated for a single hexagon in figure 2. This sign is 
encoded in (ia • cr), where 



V H = 



d 



djk x d k 3 

\ d ik x d kj\ 



(3) 



... is the three-dimensional vector (with vanishing z 
component) connecting sites % and k, and k is the 
intermediate lattice site involved in the hopping process 
from site i to site j. The vector er is defined by the Pauli 
matrices as cr = (a x , a v 1 a z ). Close to the Dirac points, the 
spin-orbit term can be written as H so ~ X so ip^ a z T z s z ij), 
where the Pauli matrices refer to the spin, sublattice, and 
Dirac point, respectively [10]. 

The last term in (1) is the Rashba spin-orbit interaction 
of strength A R , defined in terms of the spin vector s = |er, 
and the unit vector = d i -/\d i -\. The Rashba coupling, 
which is purely off-diagonal in spin (hence, spin is no 
longer conserved), breaks the z M> — z inversion symmetry, 
and arises for example, in the presence of a substrate [10]. 

For A so ^ 0, the half-filled KM model describes a 
quantum spin Hall insulator as long as A R < 2\/3A go , see 
section 3. At A R = 2\/3A go , the system is gapless, and 
increasing A R further leads to a trivial insulator [9,62]. On 
a geometry with open edges, the KM model exhibits helical 
edge states with a symmetry-protected crossing at k = 
(k = tt) for armchair (zigzag) edges. 

The KM-Hubbard (KMH) model corresponds to H KM 
with an additional, local Hubbard repulsion between 
electrons, and is written as 



2.2. Haldane-Hubbard model 

Although it breaks time-reversal symmetry and hence does 
not fall into the class of quantum spin Hall models, 
the Haldane-Hubbard model (the name is somewhat 
misleading, because the interaction is nonlocal) is included 
here for two reasons. First, the noninteracting version [65] 
has played a key role for the development of the field of 
quantum spin Hall insulators as it underlies the KM model 
[9, 10]. Second, the interacting model has been studied in 
the context of topological quantum phase transitions [66- 
68]. 

The Hamiltonian can be written as 



H 



HH 



^E 

(ii) 



nn - 

i- 3 



* 2 E eifec h 



(5) 



1) 



(4) 



The first term describes the nearest-neighbour hopping of 
spinless fermions on the honeycomb lattice. The second 
term is a spin-orbit coupling with a phase (f>.j — ±4> for 
clockwise (anticlockwise) next-nearest-neighbour hopping, 
equivalent to the spin-orbit term of the KM model ( 1 ) for a 
fixed spin direction. Also included is a repulsion V between 
spinless fermions at neighbouring lattice sites. Because 
of the restriction to spinless fermions (or, equivalently, to 
one spin sector), the model breaks time-reversal invariance, 
and for V — realizes the integer quantum anomalous 
Hall effect for any <f> at half filling [65]. Large enough 
interactions cause a quantum phase transition [66-68] to a 
charge-density-wave state, see section 4. 



It was first considered by Rachel and Le Hur [63]. 
Compared to the Hubbard model (A go = 0), the spin-orbit 
terms reduce the rotational symmetry from the hexagonal 
C 6 group to C 3 , and the spin rotation symmetry from 
5(7(2) to either 17(1) (for A R = 0) or Z 2 (for A R ^ 0). 
A detailed review of previous work on the KMH model is 
provided in sections 4 and 5, and the noninteracting case is 
considered in section 3. 

The KMH model has played an important role because 
it can be investigated with exact quantum Monte Carlo 
methods. Moreover, there is a close connection to the 
Hubbard model on the honeycomb lattice, which after 
reports of the existence of an intermediate spin liquid phase 
[57] has been at the focus of intense investigations. The 



2.3. Bernevig-Hughes-Zhang-Hubbard model 

The low-energy, continuum model for HgTe quantum wells 
involves two s and two p orbitals, each forming a Kramers 
pair, with a gap at the V point [17]. From this model, 
Bernevig et al [17] derived a two-dimensional tight-binding 
model on the square lattice that can be studied numerically. 
A possible real-space representation, hereafter referred to 
as the BHZ model, is [7] 



-^bhz = E E E e « c l 



(6) 



EEE' 
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where a labels the four bonds to nearest-neighbour lattice 
sites, a, j3 G {s, p} are orbital indices, and a is a spin index. 
The 2x2 hopping matrix defined in orbital space is given 
by (here sgn(er) = +1, -1 for a =t, |) 



t a — 



sp 



(7) 



*±v I ±isgn( ( 7)i sp 



±isgn(cr)t, 

-i. 



sin 2 fc , 



(8) 



1/2 



The momentum-space representation of the tight-binding 
model reduces to the continuum model for HgTe upon 
expanding around the V point [17]. A convenient choice of 
parameters is to set e s = — e p = e/2, so that e is the energy 
splitting between s and p orbitals, and t m — t pp = t. The 
energy bands are then given by 

E(k)= ±2[u%{s\n 2 k x 

+ [e/2 - 2t (cos k x + cos k y )] 2 1 

Writing the BHZ model in terms of Dirac matrices [7] and 
considering the eigenvalues of the parity operator at the 
time-reversal invariant points of the Brillouin zone gives a 
quantum spin Hall phase for — 8t < e < 8t, and a trivial 
insulator else. This condition, when translated into the 
corresponding notations, agrees with previous work [7, 17, 
69,70]. However, as pointed out in [69], there are in fact 
two different topological phases, separated by a metallic 
point at e = 0. These phases have a topological band gap 
at the T point, A(T) = 4|e/2 - 4t|, and at the M points, 
A(M) = 2 1 e | , respectively. They can be distinguished 
by their response to topological defects in the form of 
dislocations [69]. For systems with edges, the protected 
crossing point of the helical edge states is located at k = 
(T phase) and k — ir (M phase), respectively. Only the T 
phase is observable in HgTe quantum wells [17], whereas 
the M phase is an interesting feature of the tight-binding 
model. 

The Hamiltonian (6) neglects terms that couple the two 
spin sectors. Such contributions arise from the breaking of 
inversion and/or axial symmetry in HgTe quantum wells, 
but are expected to be small [17]. Nevertheless, the spin 
symmetric case is not generic. The effects of nonconserving 
terms in this model have been discussed, for example, by 
Schmidt et al [71]. 

Because there are two orbitals per site, the definition 
of a local Hubbard interaction term is not unique, 
and both intra-orbital interactions of the form H u — 
U 2~2ia ^io^iai [72-74] as well as intra- and interorbital 
interactions, H v = U'^2 i n i {n i — 1) (with n i = 

En- ) [70] have been studied. For realistic 
calculations, the full Coulomb interaction (including Hunds 
rule coupling) should be taken into account. Interaction 
effects on the helical edge states of the BHZ-Hubbard 
(BHZH) model have been explored in [71, 73, 74], whereas 
bulk properties were investigated in [70, 72, 73, 75]. 




it2C x it2(T y it2(T z 

Figure 3. Spin-orbit terms of the SI model (9), see also [22]. 

2.4. Sodium iridate model 

Shitade et al [22] predicted a quantum spin Hall phase with 
substantial spin-orbit and electron-electron interactions 
in the transition-metal oxide Na 2 Ir0 3 . Exploiting the 
layered structure of this material, they derived an effective, 
two-dimensional model on the honeycomb lattice. The 
Hamiltonian takes the form [22, 25, 76] 



si 



(ij) a 

^E W 



EE 4m 

((ij)) ss ' 



(9) 



It describes nearest-neighbour and next-nearest-neighbour 
electron hopping between Ir atoms. The indices s, s' refer 
to pseudospin states [22], but will also be called spin states 
in the following. The nearest-neighbour hopping (i) is real 
and diagonal in spin. The next-nearest-neighbour hopping 
in general has real and complex contributions [22], 



(10) 



where w — x,y,z selects one of the Pauli matrices 
depending on the hopping path, as illustrated in figure 3. 
The hopping t 1 is sometimes neglected. 

In contrast to the KM model (1), there is only one 
parameter t 2 that determines the strength of both spin- 
conserving and nonconserving hopping processes. Spin 
is therefore generically not conserved in the SI model. 
Moreover, both contributions act between next-nearest- 
neighbour lattice sites, whereas the spin-orbit and Rashba 
terms in the KM model (1) connect next-nearest and nearest 
neighbour sites, respectively. 

The model (9) describes a single honeycomb layer of 
Ir atoms in Na 2 Ir0 3 , and may be realized at the surface 
of a three-dimensional sample [22]. In contrast to the 
KMH model for graphene and the BHZH model for HgTe, 
strong electron-electron interactions have been argued to be 
generically present in this material. The SI model has been 
theoretically studied in [25,26,76]. Since it is defined on 
the same lattice as the KMH model, the helical edge states 
again cross at k = (k = it) on armchair (zigzag) edges. 

2.5. Other models 

Whereas a minimum of four states per unit cell is necessary 
for a gapped, time-reversal invariant band structure, there 
exist several examples of more complex models that support 
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a quantum spin Hall state. These include the kagome lattice 
[77], the decorated honeycomb lattice [78], the Lieb lattice 
[79], the square-octagon lattice [80], and the ruby lattice 
[81]. The corresponding multi-band Hamiltonians are of 
particular interest for two reasons. First, the decorated 
honeycomb lattice has been shown to provide a remarkable 
adiabatic connection between a noninteracting topological 
state (the quantum spin Hall state) and a strongly interacting 
topological state (the chiral spin liquid phase of the Kitaev 
model [82,83]), both of which are exactly solvable [35, 
78]. Second, such models can have band structures that 
include very flat bands, and thereby possess the potential 
for fractional topological states arising from electrons in 
strongly correlated, partially filled flat bands [81], see also 
section 4.5. Finally, whereas a quantum spin Hall state 
on the honeycomb lattice requires complex next-nearest- 
neighbour hopping, nearest-neighbour hopping is sufficient 
on more complicated lattices such as the kagome or the 
decorated honeycomb lattices [78]. 

Other models that are relevant in the context of 
correlation effects are due to Araujo et al [84] (a bilayer 
square-lattice model), Goryo and Maeda [85] (the bilayer 
KMH model), and Cocks et al [86] (a time-reversal 
invariant version of the Hofstadter problem). Finally, 
lattice models for topological Mott insulators and fractional 
topological insulators will be discussed in sections 4.4 
and 4.5, respectively. 

3. Noninteracting quantum spin Hall insulator 

Although this review is mainly concerned with correlation 
effects, a brief summary of the noninteracting case 
provides the starting point for the discussion of interacting 
systems. It also illustrates the key concepts of time- 
reversal invariance, the bulk-boundary correspondence, and 
topological protection. For concreteness, the example of 
the KM model (1) is used. More general reviews have been 
given in [12, 36]. 

The KM model (1) can be solved exactly by Fourier 
transformation. In the absence of a Rashba term, it can be 
written in the form H KM = Y^ k ^k H ( k )^k> witn 
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, 6jy . Here a and b operators refer 



to the two sublattices of the honeycomb lattice, see figure 1 ; 
9k = * Si e lk ' Si is related to the nearest-neighbour hopping 
(see (2)), whereas the spin-orbit coupling enters via 7 fc = 
2A so [2cos(3fc x /2) sin(V3fc a /2) - sin( v / 3fc y )]. From (11) 
it is quite obvious that the KM model can be regarded as 
two decoupled models for the t and I spins, each equivalent 



to a spinless Haldane model (5) [8] described by the 2 x 2 
matrices H rT (k). The two spin sectors become coupled in 
the presence of Rashba terms, see (1). 

Because time-reversal symmetry is inherently linked 
to the quantum spin Hall insulator, it is worth considering 
the corresponding transformation of (11). For spin-1/2 
fermions, the time-reversal operator can be written as 6 = 
cxp(iTra v )K [7,36], where a v is a Pauli matrix and K 
denotes complex conjugation. Application of 9 to a single- 
particle Bloch state leads to k >->• k and flipping of the 
spin. Interchanging the j" and 4- sectors of (11), taking the 
complex conjugate, and using j_ k = — 7 fe and g_ k = g k , 
one can verify that the KM Hamiltonian indeed has the 
symmetry H(—k) = Oi7(fe)9 _1 . In contrast, this is not 
the case for the spinless Haldane model [8]. 

The eigenvalues of (1 1) are 



E ± {k) =± 



9 k V+ll- 



(12) 



so that the spectrum has two bands, each of which has 
a Kramers degeneracy between a =t and a =|. For 
A go = 0, the band structure reduces to (2), and the familiar 
gapless Dirac spectrum is recovered at the points K, K' 
which determine the low-energy physics of the KM model. 
Equation (12) shows that any nonzero spin-orbit coupling 
A go opens a gap A go = 3\/3A gC) at the Dirac points 
[10,63]. For X SQ /t > 1/(3^3) ~ 0.2, the minimal gap 
of size A = 2t is instead located at the M points [63]. 

In the Haldane model, the spin-orbit term gives 
rise to an integer Chern invariant, a chiral edge mode, 
and a Hall conductivity a = ±e 2 /h [8]. The key 
observation of Kane and Mele [9, 10] was to realize that 
combining two copies of the Haldane model, with Hall 
conductivities = — a^. y (the opposite sign is apparent 
from the diagonal matrix elements in (11)), leads to a 
model of spinful fermions that preserves time-reversal 
symmetry. Remarkably, this construction — dubbed the 
Haldane-Kane-Mele correspondence in [35] — extends to 
interacting systems with and without intrinsic spin-orbit 
coupling, and even to fractional states [32]. Whereas the 



Hall conductivity a x 



0, as required by 



the fact that cr is odd under time reversal, the spin Hall 

x y 

conductivity a s xy = (fi/2e)(a\jj y — ^iy)' re l atec l to the spin 
current J s = (h/2e)(J^ — takes on nonzero values 
a xy = ±e/2ir [10]. The sign depends on the sign of 
the spin-orbit coupling A go [62]. The quantization of a xy 
holds as long as spin is conserved, but a nonzero spin Hall 
conductivity exists under more general conditions [10]. 

Remarkably, the quantum spin Hall phase of the KM 
model, which does not break any symmetries, represents 
a new state of matter that is topologically distinct from 
a simple band insulator [9, 10]. In close analogy to the 
integer quantum Hall effect [3, 87], a topological invariant 
can be defined. An explicit mathematical definition is most 
easily achieved for the case where spin is conserved and the 
quantum spin Hall state decouples into two integer quantum 
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Hall states, as described by (11). For a noninteracting, 
integer quantum Hall system, the Chern invariant for the 
m-th band is given by [3-5, 88] 



2tt/ 



2tt 

where F(k) = 



(13) 

A m (k) is the so-called Berry 
flux or Berry curvature (a gauge-invariant, local quantity), 
and A m — i(u m \ V fc \u m ) is the Berry phase or Berry 
connection (whose line-integral around a closed loop in 
momentum space gives the phase picked up by the Bloch 
state \u m )) [88]. A particularly simple representation can 
be obtained for a 2 x 2 Hamiltonian, such as H" in (11). 
Writing the Hamiltonian as H a (k) = h a (k) ■ er, and 
defining the unit vector field h a (k) = h a {k)/\h a (fc)|, the 
Chern number reads [89, 90] 



C a = 



An j 



^ , d 2 k d k h°{k)xd k h°{k) -h°(k). (14) 

In this case, the topological invariant corresponds to 
the winding number of the mapping h?(k) between the 
Brillouin zone and the unit sphere; mathematically, nonzero 
values of C a require three nonzero components of the 
vector h l7 (k), which physically corresponds to broken 
time-reversal invariance [90]. 

For a quantum spin Hall insulator with time reversal 
symmetry, as described by (11) for A go ^ 0, the Chern 
indices have opposite sign, so that + = 0. However, 
the spin Chern number [62,91] C s = (C* - C l )/2 = ±1 
is nonzero, and the Z 2 topological invariant defined by [36] 

i/ = C* s mod 2 (15) 

takes on the values 0,1. If spin is conserved (A R = 0), v 
is directly related to the quantized spin Hall conductivity, 
a* — ±ve/2n. Hence v = 1 corresponds to the 
topologically nontrivial state that exhibits the quantum spin 
Hall effect. When spin is not conserved, it is no longer 
possible to calculate v from the Chern indices and C^, 
and o* is no longer quantized to ±e/27r. Nevertheless, 
a quantum spin Hall phase with quantized v = 1 persists 
[9,62]. The calculation of the topological invariant of 
systems without spin conservation, and with interactions 
and/or disorder will be the topic of section 6. 

The existence of a topological state in the KM model 
relies on the special form of the spin-orbit term. The 
topological band structure arises from the opposite sign of 
the spin-orbit gap at K and K' in the presence of inversion 
symmetry. Explicitly, 7 K , = — 7 K , as can be verified using 
the definition of 7 given after (11) and the definition of 
the Dirac points in figure 1. As a counter example, a gap 
in the spectrum can also be opened by adding a staggered 
sublattice potential, X v J2 i ^ i c\c i , with = ±1 depending 
on the sublattice. Such a term, which preserves time- 
reversal symmetry but breaks inversion symmetry [8,92], 
has the same sign at K and K', and leads to a topologically 
trivial, semiconducting phase with v = [8-10]. The 
topological phase can be destroyed when the gap is closed 



by competing terms such as Rashba spin-orbit coupling [9, 
62]. 

Similar to the quantum Hall effect, the fact that the 
quantum spin Hall phase is topologically distinct from a 
simple band insulator has important consequences. As 
implied by the term topological, such as state cannot 
be connected to a nontopological state via adiabatic 
transformations of the Hamiltonian. Instead, a transition 
from a quantum spin Hall state to a topologically trivial 
state can occur either via a closing of the bulk band gap, 
or via the breaking of time-reversal invariance. 

Because a closing of the gap is required, an interface 
between a quantum spin Hall insulator and a trivial insulator 
such as the vacuum necessarily gives rise to metallic states 
[36]. The existence of these states may also be inferred 
from a flux insertion argument [93]. Whereas all electrons 
at a given edge of a quantum Hall sample move in the 
same direction {chiral edge state) [94], quantum spin Hall 
insulators have helical edge states [54]: at each edge, 
electrons with spin f move in one direction and electrons 
with spin I move in the opposite direction. Hence, there is 
a one-to-one correspondence between spin and direction of 
motion in helical edge states [10, 54]. The existence of pairs 
of edge states with opposite chirality follows directly from 
equation (11), which combines two integer quantum Hall 
states that possess chiral edge modes of opposite velocity. 

For noninteracting systems, the number of pairs of 
edge states (modulo 2) is directly linked to the value of 
the topological invariant v [9,36]. Moreover, the so-called 
bulk-boundary correspondence expresses the fact that the 
existence of edge states is guaranteed by the topological 
nature of the bulk system. This correspondence goes 
back to the work of Halperin [6]. It also follows from 
topological field theory when imposing gauge invariance 
[95], see below. Intuitively, it can be understood as a 
consequence of the protection of the topological state with 
respect to adiabatic deformations, implying that the only 
way to achieve a change of the topological invariant is by 
closing the bulk band gap. Such a closing, which can either 
occur as a function of a parameter (phase transition), or at 
the interface with another material (including the vacuum), 
is tantamount to the existence of gapless excitations. Hence, 
edge properties can in principle be studied to infer the 
properties of the bulk, but as discussed in sections 5 and 6, 
this equivalence can break down in interacting systems. 

The quantum spin Hall state cannot be destroyed 
by small perturbations such as disorder or electron- 
electron interactions [9,10,54,55]. Exploiting the bulk- 
boundary correspondence, this stability can be understood 
by considering the helical edge states. Denoting the 
band dispersion of the edge states by e a {k) (where k is 
the conserved component of the two-dimensional crystal 
momentum along the edge) time-reversal symmetry leads 
to the condition 



e 4 (-fc) = e t (fc). 



(16) 



2n;/a 




0.05 0.1 0.15 0.2 
Rashba coupling Vr ( t ) 

Figure 4. (a) Edge states of the KM model on a zigzag ribbon with open 
boundary conditions in the y direction (see inset) and AsoA = 0-03 [10]. 
(Reprinted with permission from [10]. Copyright 2005 by the American 
Physical Society), (b) Phase diagram of the KM model as a function of 
Rashba coupling Vr. (= Ar) and disorder strength W, at fixed spin-orbit 
coupling Ago/i = 0.05 [62]. The quantum spin Hall state is stable 
up to Ar, = 2\/3Aso ~ 0.17t. (Adapted with permission from [62]. 
Copyright 2006 by the American Physical Society). 



At the special, time-reversal invariant momenta k = 0, it 
(these points are mapped onto themselves under the action 
of O), this condition implies 

e t (k) = e l (k), k = 0,n, (17) 

and hence a symmetry-protected, gapless crossing point of 
edge states. It turns out that in the topological phase of 
the KM model, there exists a single pair of edge states 
that cross either at k = for an armchair edge, or at 
k = 7r for a zigzag edge [9, 10]. The latter case is 
illustrated in figure 4(a) for a ribbon with open boundaries 
in the y direction, whereas the armchair case can be seen in 
figure 16. The edge states in figure 4(a) connect the Dirac 
points K and K', and their helicity becomes apparent by 
plotting the spin-resolved single-particle spectrum. 

Because time-reversal symmetry implies that the edge 
state dispersions are degenerate at k = or it [10], pertur- 
bations that preserve time-reversal symmetry cannot open 
a gap. Most importantly, single-particle backscattering, 
which is the most relevant perturbation in ordinary one- 
dimensional metals [96], is not allowed in a helical liq- 
uid because the single-particle states |fc,t) and |— k,\.) are 



Kramers partners, and therefore orthogonal. Consequently, 
provided time-reversal symmetry is not broken, the heli- 
cal edge states are protected from elastic, single-particle 
backscattering [10, 12,54,55], but the protection is weaker 
than that of chiral edge states in a quantum Hall insulator 
where all electrons at a given edge propagate in the same di- 
rection, and there is no phase space for backscattering. The 
presence of counter-propagating electrons at a given edge 
has important implications for strongly interacting systems, 
see section 5. Finally, an even number of edge states exists 
for each spin direction in the trivial phase (v = 0) of the 
KM model, and scattering is not forbidden by symmetry 
[54,55]. 

The stability of the edge states with respect to disorder 
and Rashba spin-orbit coupling has been numerically 
verified [9,62,97], and the resulting phase diagram is 
reproduced in figure 4(b). The topological phase exists 
as long as the Rashba coupling does not exceed A R = 
2v / 3A go pa 0.17t (at this point, the topological gap closes 
and the system becomes a trivial insulator), and survives 
up to a critical value of the disorder strength. This result 
is crucial for the experimental realization of the quantum 
spin Hall effect [14, 15] because neither Rashba coupling 
nor disorder can be completely excluded. In the presence 
of interactions, the stability of the edge states is no longer 
guaranteed, as discussed in section 5. 

As shown in this section, many aspects of integer 
quantum Hall and spin Hall insulators can be understood in 
the framework of topological band theory [12]. However, 
the latter is inherently restricted to noninteracting systems. 
An alternative approach to topological states of matter is 
based on Chern-Simons field theory, a low-energy theory 
of gauge fields that is applicable to interacting and even 
fractional states. A comprehensive introduction to Chern- 
Simons theory can be found, for example, in [2]. For 
the extension to topological insulators with time-reversal 
symmetry see [53]. 

To illustrate how Abelian Chern-Simons theory 
explains the quantization of a xy as well as the presence 
of edge states in a quantum Hall insulator, it is sufficient 
to consider the noninteracting Haldane model in the form 
H = Y^tij *i-j c ] c i> c f- (5), and couple it minimally to an 



electromagnetic field, 

H(A") = XX-cj^.exp 



+ ec 



27ri 



A(l,t) -dl 



(18) 



The scalar and vector potentials may be combined in the 
form = (A Q , -A) with A = $/c and A = (A 1 ,A 2 ). 
Gauge invariance implies — with j M = —SH/SA^. 
Linear response theory then directly gives the quantized 
Hall conductivity. 

The starting point for the derivation of Chern-Simons 
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theory is the Grassmann path integral [34] 

J i 

with the action 



S(A U ) 



at 



<*(«) 



(19) 



(20) 



f(i,t)A(i,t) 



Equation (20) follows from a gradient expansion in A^. The 
fermions can be integrated out, at the expense of a Gaussian 
integral, to obtain Z(A ) = e iScff ( A »^ with 



(21) 



Here P^ v is the current-current correlation function, which 
is short-ranged in an insulating state, and a summation 
over repeated Greek indices is implied. Hence, to a first 
approximation, retardation effects can be omitted by setting 
t' = t. Taking the continuum limit, all possible gauge 
invariant terms can be derived. The relevant contribution 
for topological states is the Chern-Simons term 

„2 



4-7T 



(22) 



from which all the low-energy properties of the integer 
quantum Hall state follow. The invariance of the theory 
with respect to twists of the boundary conditions by a 
multiple of the flux quantum <i> pins the first Chern 
number C to an integer [3,4]. Using j = 5S^ /5A^, 
and assuming A Q = 0, gives the familiar quantized Hall 
response 

„2 



3: 



= a xy E y > 



■>■!) 



= c- 



(23) 



On a torus geometry, the Chern-Simons action (22) is 
invariant under the gauge transformation 4 — > 4 + 9 v. 
However, in the presence of a boundary (e.g., on a cylinder), 
the gauge transformation leads to an additional boundary 
term that is directly related to the chiral edge state of a 
quantum Hall insulator [2]. 

Finally, a Chem-Simons theory of a noninteracting, 
spin-conserving quantum spin Hall insulator (as appropriate 
for the KM model) takes the form 

„2 



cCS _ 
J QSH — 



47T 



dt 



For more general cases, including fractional topological 
states, see [98] and references therein. 



i xyz AFMI 



xy AFMI 
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Figure 5. Phase diagram of the half-tilled KMH model from quantum 
Monte Carlo simulations [99-101]. The phases correspond to a semimetal 
(SM, for Ago = 0), a quantum spin Hall insulator (QSHI), a quantum 
spin liquid (QSL), and an antiferromagnetic Mott insulator (AFMI) with 
either Heisenberg (for Ago = 0) or easy plane (for Ago 7^ 0) order. Here 
A R = 0. Data taken from [99-101]. 



4. Bulk correlation effects 

The integer quantum spin Hall effect can be understood by 
considering only the noninteracting band structure. Much 
of the recent theoretical work has focused on the impact 
of electronic correlations. Given the scope of this review, 
the discussion of bulk correlation effects begins with results 
for models with intrinsic spin-orbit coupling that support 
an integer quantum spin Hall state in the noninteracting 
limit. The investigation of such systems is particularly 
rewarding as the corresponding models can be studied using 
numerical methods, including exact quantum Monte Carlo 
simulations that permit to investigate, for example, the 
critical behaviour at quantum phase transitions. After an 
overview of phenomena using the example of the intensely 
studied KMH model, a detailed discussion of correlated 
quantum spin Hall phases and interaction-driven phase 
transitions with or without symmetry breaking is given, 
followed by two shorter sections on interaction-driven 
topological insulators, and fractional topological insulators. 

4.1. Phases of the Kane-Mele-Hubbard model 

The KMH model provides a framework to study the effects 
of electronic correlations on a quantum spin Hall insulator. 
It is of particular interest due to connections with graphene 
[59] and the Hubbard model on the honeycomb lattice 
[57]. Given its historical importance and the possibility of 
applying exact quantum Monte Carlo methods, the KMH 
model is arguably the best understood model among those 
introduced in section 2. 

Figure 5 shows the numerically obtained [99-101] 
phase diagram of the KMH model at half filling, and 
with conserved spin (A R = 0). Most importantly, the 
numerical results confirm the stability of the topological 
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state of the KM model (U = 0) up to Hubbard interactions 
U/t > 4, as first predicted based on mean-field theory [63]. 
As discussed in more detail below, the numerical results 
suggest that the topological phase at U > is adiabatically 
connected to the noninteracting case. It is therefore referred 
to as a topological band insulator. Even stronger electron- 
electron repulsion drives a quantum phase transition to a 
magnetic state discussed in detail in section 4.3. The 
other phases of the KMH model shown in figure 5 are 
a semimetal (SM) that exists for vanishing spin-orbit 
coupling (with a massless Dirac spectrum similar to (2)), 
and possibly a gapped quantum spin liquid (QSL) phase 
at intermediate U jt [57] to be discussed in section 4.3. 
Contrary to analytical predictions [102], numerical studies 
so far indicate that all phase transitions in figure 5 are 
continuous [100]. 

The phase diagram of the KMH model has been 
investigated by a variety of theoretical methods. Mean-field 
theory captures the magnetic transition, and the increase of 
the critical U with increasing A go [63]. Quantum Monte 
Carlo simulations provide the exact phase boundary for 
the magnetic phase, demonstrate the absence of magnetic 
order in the z direction, and suggest the existence of a 
quantum spin liquid phase [99, 100, 103]. (Figure 17 shows 
the numerical phase diagram obtained by Zheng et al 
[103].) Approximate variational Monte Carlo results are 
also available at A so /i = 0.1 [104]. Quantum cluster 
methods such as the cluster dynamical mean-field theory 
[105] and the variational cluster approach [106, 107] are 
able to capture the overall structure of the phase diagram. 
However, care has to be taken regarding the choice of 
the cluster shape and size. Analytical approaches include 
slave-boson [108], slave-spin [109] and slave-rotor methods 
[70], as well as mappings to spin models to explore the 
regime of large U/t [25,26, 110]. Using field theory and 
symmetry arguments, predictions about the phases and 
phase transitions of the KMH model were made by Griset 
and Xu [102] as well as by Lee [111]. 

With the exception of the quantum spin liquid phase, 
whose existence is still under debate [57,99, 100, 112] and 
expected to be specific to the honeycomb lattice, the phase 
diagram in figure 5 contains many features that are generic 
for models of correlated quantum spin Hall insulators. The 
semimetal, although not a topological state, also exists in 
the SI model defined in (9). In fact, in the absence of 
next-nearest-neighbour hopping (A go = in (4) and t x = 
tn = in (9)), both models reduce to the Hubbard model 
on the honeycomb lattice. A quantum phase transition to 
a magnetic state with easy-plane order (as a result of spin- 
orbit coupling) at large U/t has also been observed in the 
BHZH model [70,72], the SI model [76], and a spinful 
version of the Hofstadter butterfly problem with Hubbard 
repulsion [86]. The existence of a magnetic phase is a direct 
consequence of the nonfrustrated exchange interactions that 
emerge on bipartite lattices. Similarly, a correlation-driven 



0.6 



,0.4 



O A sp , U/f=2(QMC) 
□ A ,U/f=2(QMC) 
_ A,l//f=0 

sp 

-- A ,U/t=0 




0.02 0.03 0.04 0.05 
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Figure 6. Single-particle gap A ap and spin gap A s as a function of 
spin-orbit coupling A = Ago from quantum Monte Carlo simulations 
of the KMH model with U/t = 2 [100]. Also shown are the gaps of 
the noninteracting KM model (see text). (Reprinted with permission from 
[100]. Copyright 2012 by the American Physical Society). 



transition from a topological to a magnetic state is also 
expected in three dimensions [29]. 

4.2. Correlated topological band insulators 

In section 3, the bulk-boundary correspondence was 
invoked to argue that the quantum spin Hall phase should be 
stable as long as time-reversal invariance is not broken. This 
idea is borne out by the phase diagram in figure 5, which 
shows that (with the possible exception of the region around 
the spin liquid phase) the topological state is eventually 
destroyed when time-reversal symmetry is spontaneously 
broken. Such behaviour seems to be typical of models that 
conserve spin [25]. On the other hand, an interaction-driven 
transition to a fractional topological state that conserves 
time-reversal symmetry has been reported for the SI model 
[76] (see section 4.3). 

Given a correlated quantum spin Hall insulator that 
exists up to rather large values of U/t, an important 
question is to what extent this phase resembles that of the 
noninteracting KM model. A partial answer is provided by 
the quantum Monte Carlo results for the single-particle gap 
A sp and the spin gap A g shown in figure 6. Starting in 
the correlated semimetal phase at U/t = 2, any nonzero 
spin-orbit coupling A go opens both gaps simultaneously. 
Importantly, for small A so , both gaps follow very closely 
the corresponding results for the U — KM model 



(indicated by the lines in figure 6), namely A sp = A 



so 



3y3A and A s = 2A sp . This observation suggests that the 
quantum spin Hall states with U/t — 2 and U/t = are 
very similar. Minor correlation effects are apparent from 
the fact that the gaps for U/t = 2 are systematically smaller 
than those for U/t = 0, and that A s < 2A sp . 

Yoshida et al [72] have used dynamical mean-field 
theory to calculate the spin Hall conductivity of the BHZH 
model as a function of temperature and the intra-orbital 
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Figure 7. Spin Hall conductivity as a function of interaction strength 
and for different temperatures, as obtained from dynamical mean-field 
calculations for the BHZH model with intra-orbital repulsion U [72]. 
(Reprinted with permission from [72]. Copyright 2012 by the American 
Physical Society). 

repulsion U. The results, shown in figure 7, reveal a 
clear suppression of a s xy due to thermal fluctuations and 
electronic correlations, and that at,, eventually vanishes 
at U/t ~ 12. However, as demonstrated in the inset, 
for small enough U /t, the conductivity extrapolates to the 
quantized value e 2 /h for T —> 0. Consequently, although 
interaction effects are clearly visible at finite temperatures, 
the correlated quantum spin Hall phase of the BHZH 
model resembles the noninteracting state at T = 0. As 
pointed out in section 3, a nonquantized Hall conductivity 
is also observed in the noninteracting case when spin is 
not conserved; however, the topological invariant retains 
its integer value v — 1 [62]. It would be interesting to 
confirm this explicitly in an interacting model, see also 
section 6. Temperature effects on the quantum spin Hall 
phase of the KMH model have been discussed in [105]. 
The topological nature of the quantum spin Hall phase at 
U > has also been demonstrated for the KMH model 
both at T = and T > with the help of ir fluxes [101]. 
Finally, the topological invariant of the KMH model has 
been calculated using the variational cluster approach [107]. 

The properties of the correlated quantum spin Hall 
phase have further been explored by numerically calcu- 
lating the single-particle spectrum in systems with edges, 
see section 5. Both for the KMH model [105, 106] and 
the BHZH model [73], gapless helical edge states are ob- 
served even for substantial values of U/t, although interac- 
tion effects modify the spectral weight, edge transport and 
magnetic correlations (see section 5). In summary, in the 
spin-conserving KMH model, the quantum spin Hall phase 
remains qualitatively unaffected by correlations up to the 
point where the state breaks down, suggesting an adiabatic 
connection to the noninteracting spin Hall insulator. This 
finding has been exploited to construct an effective model 
of the edge states in which a Hubbard interaction is only 
taken into account at the edge, whereas the noninteracting 
bulk establishes the topological phase [99, 113]. 



4.3. Interaction-driven phase transitions 

Although the quantum spin Hall phase of, for example, 
the KMH model persists even in the presence of rather 
strong interactions, it eventually breaks down at the onset of 
long-range magnetic order that spontaneously breaks both 
the U(l) spin symmetry and time-reversal symmetry [63, 
99, 103]. Such a magnetic transition is expected to be a 
generic feature of quantum spin Hall models. In the KMH 
model, it can be studied in detail by quantum Monte Carlo 
simulations. Whereas the magnetic transition involves a 
well-defined order parameter and can be captured already 
in mean-field theory, interactions can also lead to more 
complex phenomena in which quantum fluctuations play a 
key role. These include topological Mott insulators [28], 
quantum spin Hall states with Mott insulating edges [29, 
63], or quantum spin liquids [99]. It is even possible to 
have a coexistence of a topological insulator with protected 
edge states and long-range magnetic order [1 14]. 

4.3.1. Symmetry-breaking transitions The best studied 
example of an interaction-driven quantum phase transition 
that involves the breaking of time-reversal symmetry is the 
magnetic transition observed in the KMH model. As visible 
from the phase diagram in figure 5, the Hubbard repulsion 
induces a transition from the quantum spin Hall phase to a 
Mott insulator with long-range magnetic order at a critical 
value Ujt. The existence of this magnetic transition has 
first been demonstrated using mean-field theory [63]. 

Already for U < U , inside the correlated quantum 
spin Hall phase, the Hubbard interaction gives rise to 
the formation of local moments and substantial magnetic 
correlations. These correlations can be measured by the 
transverse spin structure factor [99, 100, 103] 

OL 

ra^ = p£(-i) a (-iy' 

rr' 

and, similarly, by the longitudinal spin structure factor 

sir = Ej^r with 

[s^r = ^£(-i) a (-i) <*oi^ a ^i*o} • ( 26 ) 

rr' 

Here r, r' are vectors that indicate unit cells of the 
honeycomb lattice, a, /3 € {A, B} are sublattice indices, 
(-1)" = 1 (-1) for a = A (B), and the trace of the 2x2 
matrices is taken. 

Quantum Monte Carlo results for the quantity ST^/N 
(with N = 2L 2 being the number of lattice sites), which is 
related to the transverse magnetization by m xy = S x ^ v /N, 
are shown in figure 8 for different values of U/t and 
for X SQ /t = 0.1 [100]. Up to U/t = 4.9, magnetic 
correlations are substantial on small length scales but the 
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Figure 8. Finite-size scaling of the rescaled magnetic structure factor 
S^ F /N (25) in the KMH model with X S o/t = 0.1 and for different 
values of U/t [100]. The quantity S^L/N extrapolates to zero in the 
quantum spin Hall phase and to a finite value in the magnetic phase; the 
critical value is U c /t = 4.95(5). Inset: order parameter m xy as obtained 
from the size-extrapolated values. (Reprinted with permission from [100]. 
Copyright 2012 by the American Physical Society). 

structure factor extrapolates to zero in the thermodynamic 
limit L — > oo. For U/t > 5, long-range magnetic 
order is found. The critical value can be estimated as 
Ujt = 4.96(4) [100], and the transition is well visible 
from the magnetization results displayed in the inset of 
figure 8. A similar analysis for SV/ ¥ /N reveals that no long- 
range order exists up to very large values of U /t [100], in 
contrast to predictions of longitudinal order based on cluster 
calculations [106]. A very similar picture arises for other 
values of X so /t [100]. At first sight, the increase of U c 
with increasing A so visible in figure 5, and also observed 
in several other works [63, 70, 99, 103, 105-107, 1 10], may 
be attributed to the increase of the spin-orbit gap. However, 
since U c also increases from Ujt — 5.70(3) at X so /t = 
0.2 [101] to Ujt = 6.5(1) at X so /t = 0.3 [99] even 
though the minimal gap at the M points remains constant 
[63], the origin of the increase of U is instead expected to 
be related to the competition between kinetic energy and 
magnetic order. 

The fact that the magnetic ordering occurs in the xy 
direction can be understood from an effective spin model 
derived in the limit U/t — > oo where charge degrees of 
freedom are frozen out [63]. The spin Hamiltonian takes 
the form [25,63, 110] 

#kmh = Jj2 s i- S i+ J ' E {s;s]-s?s?-s?S3) .(27) 

The first term is an antiferromagnetic coupling between 
neighbouring lattice sites emerging from nearest-neighbour 
hopping (J = 4t 2 /U), whereas the second term originates 
from the spin-orbit interaction (J' = 4Ag /?7) and couples 
next-nearest-neighbour sites ferromagnetically in the xy 
direction and antiferromagnetic ally in the z direction. 
Hence, at least for large U / 1, the KMH model corresponds 
to a frustrated spin model with respect to the z component 
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Figure 9. (a) Rescaled transverse magnetic structure factor S^p / ^> see 
(25), as a function of U in the KMH model with Aso/* = 0.1. Assuming 
3D XY behaviour and using the critical exponents of the 3D XY model 
[1 15], a clean intersect of curves for different system sizes L at the critical 
point U c /t = 4.96(4) is obtained, (b) Scaling collapse of L 2 ' 3 /"5^ F /iV 
as a function of L 1 / I/ ((7 - U c )/U c using U c /t = 4.96(4) and 3D XY 
exponents. (Reprinted with permission from [100]. Copyright 2012 by the 
American Physical Society). 

of spin, but allows for nonfrustrated magnetic order in the 
direction. Exact numerical results [99, 100, 103] confirm 
this strong-coupling picture. For the parameter range of 
figure 5, no order in the z direction was found down to 
A SQ /t = 0.002 [100]. 

Considering the spin-conserved case (A R = 0), the 
spin-orbit coupling in the KMH model reduces the spin 
symmetry from SU(2) to U(l). Therefore, the magnetic 
transition is expected to be in the universality class of the 
3D XY model [99, 102, 111]. The quantum Monte Carlo 
method permits to test this conjecture. Assuming that the 
transition is of the 3D XY type, the order parameter S^ F /N 
is expected to fulfil the scaling equation [100] 

S^/N = L- 2 ' 3 ^f 1 [(U - U C )L^ U ] (28) 

where z = 1, v = 0.6717(1) and j3 = 0.3486(1) are the 
critical exponents of the 3D XY model [115]. In particular, 
(28) implies that at the critical point, U = U c , S°/^/N 
should become independent of the system size L. As 
illustrated in figure 9(a) for X so /t = 0.1, the numerical 
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Figure 10. Single-particle gap A sp and spin gap A s of the KMH model 
as a function of U at A = Aso = O.lt [100]. At the transition 
(U c /t = 4.95(5)), A sp shows a dip, whereas A s closes. Inset: A sp and 
magnetization from mean-field theory. (Reprinted with permission from 
[100]. Copyright 2012 by the American Physical Society). 

results [115] indeed reveal the expected behaviour, showing 
a well-defined intersect of curves for different L at the 
critical point. Moreover, rescaling the U axis according 
to (28), a clean scaling collapse is obtained. The 3D 
XY scaling of S^/N has also been verified for other 
values of A go /i [100, 101], as well as for the spin gap 
A s [100]. The 3D XY criticality can be expected to be 
generic for quantum spin Hall models with conserved spin. 
If spin is not conserved, the symmetry is reduced to the Z 2 
Ising universality class, allowing magnetic order at finite 
temperatures. 

At the magnetic transition, time reversal symmetry — 
underlying the topological protection of the quantum 
spin Hall state — is broken. Therefore, a change of the 
topological invariant from nontrivial to trivial does not 
require the closing of any gaps. Indeed, quantum Monte 
Carlo simulations reveal only a slight dip of the single- 
particle gap near the critical point, see figure 10. Very 
similar behaviour is also found in mean-field theory, see 
inset of figure 10. On the other hand, the spin gap A s 
does close at the magnetic transition, corresponding to the 
condensation of magnetic excitons [100,101,111]. For 
U < U c , excitons have a finite energy, and can be used 
in combination with it fluxes (see section 6) to construct 
quantum spin models with tunable interactions [101]. The 
coupling between the helical edges states and the magnetic 
excitons at the critical point was recently discussed by 
Grover and Vishwanath [116]. 

An interaction-driven magnetic transition has also 
been reported for the BHZH model [70,72,73]. Whereas 
existing work has focused on the case of equal bandwidths 
for hopping via the s and p orbitals (t as = i ), a much 
richer picture can be expected for more general parameters. 
In the absence of the complex spin-orbit coupling (t = 0), 
the BHZH model reduces to the two-band Hubbard model. 
The latter has been studied intensely as it may describe 



consecutive Mott transitions within the heavier and lighter 
band [117-120]. The interplay of orbital-selective Mott 
transitions and spin-orbit coupling represents a fascinating 
topic for future research. The SI model (9) is expected 
to undergo a magnetic transition in the Ising universality 
class if the spin-orbit coupling is not too large [76], see also 
figure 11. 

Time-reversal symmetry is a key concept for the 
understanding of topological insulators, and underlies many 
of their properties. Remarkably, topological insulators can 
exist even in systems where time-reversal symmetry is 
broken by long-range magnetic order. As first discussed 
by Mong et al [114] for the three-dimensional case, a 
topological state with protected edge states can exist in 
the absence of time-reversal symmetry if S — QT X , 2 
remains conserved, where T x , 2 is a translation symmetry 
defined by the unit cell of magnetic order. A simple 
example is given by an antiferromagnet in which the 
magnetic order leads to a doubling of the unit cell. Three- 
dimensional antiferromagnetic topological insulators may 
be realized in GdBuOt [114] as well as with ultra- 
cold atoms [56]. In two dimensions, the existence of 
edge states protected against disorder has been shown to 
additionally require spin conservation [121]. According 
to Guo et al [121], nonmagnetic or magnetic staggered 
potentials can drive the noninteracting BHZ model into 
a topological phase. Surprisingly, for the same model, 
a correlated topological phase with coexisting magnetic 
order has recently been reported by Yoshida et al [75]. 
If this phase indeed exists, it is interesting to understand 
why it is absent in the KMH model. However, another 
possibility is that the coexistence is due to the fact that 
the dynamical mean-field theory applied in [75] does not 
allow for magnetic order in the xy direction (as preferred 
for systems with spin-orbit coupling). Indeed, at the level 
of static mean-field theory, two different critical points can 
be observed when decoupling the interaction in terms of the 
magnetization in the z and the xy direction, respectively. 
Phases in which topological order coexists with magnetic 
order, for example topological spin density waves, have 
also been found in the interacting, spinful Haldane model 
(also known as the topological Hubbard model) [122-124]. 
In contrast to the BHZH model of [75], the topological 
Hubbard model breaks time-reversal symmetry from the 
outset. 

A related theoretical problem, yet with important 
differences, is that of the Haldane-Hubbard model defined 
in (5). Without spin-orbit coupling (t 2 = 0), the nearest- 
neighbour repulsion drives a second-order transition to a 
charge-density-wave insulator that spontaneously breaks 
inversion symmetry [28, 125, 126]. The two degenerate 
ground states have all electrons on either the A or the 
B sublattice. Because inversion symmetry is in principle 
not relevant for the stability of the topological state of 
the Haldane model, in contrast to the role of time-reversal 
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symmetry for the quantum spin Hall state, the question 
how the transition will occur for nonzero t 2 is very 
interesting. This problem was studied using mean-field 
theory, exact diagonalization, and constrained-path Monte 
Carlo simulations [66-68]. Although the results do not 
fully agree, there is evidence that the Chern number changes 
from nontrivial to trivial, and that this change and the onset 
of long-range charge order can either coincide, or happen in 
two distinct transitions. It is not yet clear if a gap closes at 
the critical point. However, because the topological phase 
does not rely on inversion symmetry, and the topological 
protection is therefore not destroyed by the charge order, 
one may expect to see the closing of an excitation gap at the 
point where the Chern number changes. This appears to be 
confirmed by recent results [68], which relate the closing of 
the first excitation gap (E 1 — E , the difference between the 
two lowest eigenvalues) to a topologically protected level 
crossing that exists even in finite systems. 

4.3.2. Symmetry-conserving transitions The magnetic 
and charge-density-wave transitions discussed above in- 
volve a well-defined, local order parameter, and in general 
may or may not coincide with the destruction of topolog- 
ical order (in the sense of a nonzero Chern or spin Chern 
number, respectively). Here the possibility of interaction- 
driven transitions without symmetry breaking is discussed, 
including transitions across which topological order either 
remains or emerges. 

In the context of the KMH model, the possible 
quantum spin liquid phase shown in figure 5 has received 
a lot of attention. Quantum Monte Carlo results suggest 
that this phase has nonzero single-particle and spin gaps, 
but does not break any symmetries [57]. Furthermore, 
numerical results are compatible with a closing of the 
single-particle and spin gaps at the transition to the 
quantum spin Hall state, which suggests that the two 
phases are not adiabatically connected [99, 100]. However, 
despite the remarkable system sizes accessible in numerical 
simulations, the length scales for spin correlations are still 
not under control [100, 112], making extrapolations to the 
thermodynamic limit very delicate. The existence of the 
quantum spin liquid phase in the Hubbard model (A so = 
0) has recently been challenged based on quantum Monte 
Carlo results on lattices with up to 36 x 36 unit cells [112]. 
Signatures of the spin liquid phase in the KMH model 
have also been reported from quantum cluster simulations 
[105, 106], although the latter cannot describe a true spin 
liquid phase because of the inherent breaking of translation 
symmetry; instead, the observed paramagnetic insulator 
may be regarded as a valence bond solid. Furthermore, 
the quantum spin liquid phase of the KMH model has been 
related to a topological Mott insulator [108]. Because the 
KMH and SI models coincide for A go = t 2 = 0, the 
same phase may also be expected to exist in the SI model. 
At finite spin-orbit coupling, the frustrated spin exchange 



interactions of the SI model (see below) may even provide 
a more favourable setting for such a state. A review of the 
connections between spin liquids and topological insulators 
has been given by Fiete et al [35]. 

For two-dimensional models with explicit spin-orbit 
coupling, the possibility of a topological Mott insulator has 
first been discussed by Rachel and Le Hur [63], following 
an earlier suggestion by Pesin and Balents [29]. Within 
a slave-rotor approach to the KMH model, the charge 
degrees of freedom can be gapped out, whereas the spin 
degrees of freedom inherit the topological band structure, 
and hence preserve time-reversal symmetry [63]. Such an 
interaction-generated state is expected to have topological 
order. For the KMH model, the corresponding phase was 
suggested to lie in between the quantum spin Hall and 
antiferromagnetic Mott insulator regions of figure 5, but 
is not visible in numerical simulations [100]. Its absence 
can be attributed to the role of gauge fluctuations in two 
dimensions [63]. A topological Mott insulator with gapped 
charge and spin excitations in the KMH model [108] — 
presumably stable with respect to gauge fluctuations — was 
argued to be related to the quantum spin liquid phase 
shown in figure 5. Transitions from topologically trivial to 
nontrivial states upon increasing the Hubbard repulsion U 
have also been reported for the BHZH model [73] and for 
a bilayer square-lattice model [84], in both cases starting 
with band parameters outside the respective quantum spin 
Hall phase. 

A particularly interesting, strongly correlated topolog- 
ical phase has been suggested to exist in the SI model (9). 
Following the suggestion of a strongly correlated quantum 
spin Hall phase resulting from the coupling to a dynamical 
Z 2 gauge field by Ran et al [ 1 27] , Ruegg and Fiete [76] used 
a self -consistent Z 2 slave-spin mean-field theory to demon- 
strate the existence of the so-called QSH* phase in the SI 
model. Their phase diagram, reproduced in figure 11, ex- 
hibits an extended regular QSH phase, and a valence bond 
solid at large U jt. The latter is expected to become an an- 
tiferromagnetic Mott insulator (AFI) phase similar to fig- 
ure 5 when adding the possibility of time-reversal symme- 
try breaking to the theory [76]. For large spin-orbit cou- 
pling t 2 and large Hubbard U, the QSH* phase emerges. It 
has been pointed out that because the QSH-QSH* transi- 
tion corresponds to an order-disorder transition of the slave 
spins [76], the QSH* phase should be considered an orthog- 
onal metal rather than a Mott insulator [128]. 

According to [76], the QSH* phase is characterized 
by a bulk band gap, and protected helical edge states for 
charge and spin excitations. It supports fractional, bosonic 
quasi-particle excitations with semionic mutual braiding 
statistics [76]. The phase does not break any symmetries 
(and should hence be regarded as a spin liquid rather than a 
quantum spin Hall insulator), but has a four-fold topological 
degeneracy on a torus [76]. It is thus topologically ordered 
in the sense of section 1; fractional time-reversal invariant 
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Figure 11. Mean-field phase diagram of the SI model, showing a quantum 
spin Hall phase (QSH), a valence bond solid (VBS), and a fractional QSH* 
phase [76]. (Adapted with permission from [76]. Copyright 2012 by the 
American Physical Society). 

liquids will be further discussed in section 4.5. Because 
of the strong interactions, the edge states are expected to 
be strongly correlated [76], leading to a suppression of the 
Drude weight in agreement with results for other models 
(section 5.2). 

The reliability of the mean-field results in [76] is not 
yet clear. Moreover, because spin is not conserved, the SI 
model can presently not be investigated by quantum Monte 
Carlo methods. This situation has motivated two recent 
studies of the limit U/t = oo, in which an effective spin 
model can be derived [25, 26]. It has the form [25, 26] 

(ij) ((H)) 

-J 3 J2(S i -S j -2S?S™). (29) 

((ij)) 

Here the exchange constants are J 1 = 4t 2 /U, J 2 = 
At x 2 /U, and J 3 = 4t 2 /U. The coupling J 2 was neglected 
in [25]. Importantly, the term Sf S™ depends on the link 
that connects sites i and j, leading to frustration and the 
absence of any preferred direction for magnetic order [25]. 
In contrast, in the U/t = oo limit of the spin-conserving 
KMH model, (27), order in the spin-z direction is frustrated, 
but easy-plane order is not. The Hamiltonian (29) has close 
connections to the Kitaev-Heisenberg model [23,24], an 
exactly solvable model with a spin liquid ground state. 

For small values of J 2 , the model (29) undergoes 
transitions from stripe order to Neel order to spiral order 
upon increasing J 3 from negative to positive values [25, 
26]. For the repulsive SI model, the transition from Neel to 
spiral order is of particular interest, and has been interpreted 
as a strong-coupling signature of the VBS (or AFI) to QSH* 
transition [25,26]. A magnetic transition from the QSH* 
to the spiral state is expected to occur with increasing U / 1 
[25]. The spin-nonconserving SI model hence appears to 
describe phenomena beyond the numerically studied KMH 
model with A R = 0, although incommensurate order 



has also been predicted for the latter analytically [110]. 
An interesting question for future research is if fractional 
topological states and exotic magnetic states can also occur 
in the KMH model with A R ^ [25]. The corresponding 
spin Hamiltonian has been derived in [25]. 

4.4. Interaction-driven topological insulators 

Spin-orbit coupling plays a crucial role for topological 
insulators [9, 10, 17]. In general, materials with heavy 
atoms are required for the effects of spin-orbit coupling 
to be observable. On the theoretical side, the existence 
of topological phases in models with intrinsic spin- 
orbit coupling can be understood in the framework of 
topological band theory [12] without considering electronic 
correlations. A completely different approach to create 
quantum spin Hall phases is based on the theoretical 
concept of a topological Mott insulator, a state with an 
interaction-generated band gap and protected edge states 
[28]. In addition to providing a much richer theoretical 
setting to study correlated topological phases, this concept 
promises the existence of topological phases in a wider 
range of materials by abandoning the need for strong 
intrinsic spin-orbit coupling. 

Raghu et al [28] considered an extended Hubbard 
model on the honeycomb lattice with the usual nearest- 
neighbour hopping term (see (1)) as well as onsite (U), 
nearest-neighbour (V,) and next-nearest-neighbour (V 2 ) 
repulsion in the form of density-density interaction terms. 
In the spinless case, the repulsion V 2 frustrates the 
charge-density-wave order driven by V v and a mean- 
field treatment suggests that large values of V 2 lead to a 
state with spontaneously broken time-reversal symmetry 
characterized by a nonzero, imaginary value of the bond 
order parameter \ij = (c]cj) [28]. This state is fully 
equivalent to the integer quantum Hall state of the Haldane 
model [8]. In the spinful case, the nonzero bond order 
parameter can either have the same or opposite signs for 
the two spin directions. The former case leads to a quantum 
anomalous Hall state, whereas the latter case corresponds 
to a quantum spin Hall state from dynamically generated 
spin-orbit coupling that breaks the SU (2) spin symmetry 
and has protected helical edge states [28, 129]. At the mean- 
field level, the quantum Hall and quantum spin Hall states 
are degenerate. Taking into account quantum fluctuations 
within the random-phase approximation reveals that the 
quantum spin Hall state is favoured because of the existence 
of a Goldstone mode, and leads to the phase diagram 
shown in figure 12 [28]. In addition to the quantum spin 
Hall phase, the latter contains semimetal, spin-density- 
wave and charge-density-wave states. Similar results and 
the possibility of a spin-dependent Kekule order parameter 
were discussed by Weeks and Franz [130]. In contrast 
to models with intrinsic spin-orbit coupling, the quantum 
spin Hall state from dynamical spin-orbit interactions is 
characterized by a local order parameter. Nevertheless, it 
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Figure 12. Mean-field phase diagram of the extended Hubbard model 
(with repulsions U, Vi, V2) on the honeycomb lattice [28]. The phases 
correspond to a semimetal (SM), charge-density-wave order (CDW), 
spin-density-wave order (SDW), and a quantum spin Hall state (QSH) 
emerging from dynamically generated spin-orbit coupling. (Reprinted 
with permission from [28]. Copyright 2008 by the American Physical 
Society). 

will in general be necessary to calculate the topological 
invariant to distinguish such a state from one that has 
a nonzero order parameter but is topologically trivial. 
Whether the quantum spin Hall state from dynamical 
spin-orbit coupling can be adiabatically connected to a 
noninteracting QSH state (or instead supports fractional 
excitations) is an open question. 

Interaction-driven topological states in models without 
intrinsic spin-orbit coupling have also been studied on 
the checkerboard [64], kagome [131,132], and decorated 
honeycomb lattices [132]. Most importantly, these more 
complicated geometries support both Dirac band crossing 
points and quadratic band crossing points depending on the 
band filling in the noninteracting limit. For example, the 
low-energy physics on the kagome lattice is determined 
by a Dirac point for n = 1/3 and by a quadratic 
crossing point for n = 2/3. Both filling fractions support 
a quantum spin Hall phase in the presence of intrinsic 
spin-orbit coupling [77]. However, whereas topological 
order due to interaction-generated spin-orbit terms is the 
leading instability for n = 2/3, it competes with other 
ordering phenomena and requires strong interactions and 
fine-tuning in the case n = 1/3 due to the vanishing 
density of states at the Fermi level [132]. In principle, 
because these lattices support partially filled flat bands, 
the interacting models are good candidates for fractional 
topological phases. A topological phase arising from 
spontaneous ordering of complex orbitals has been found 
theoretically in a model of transition-metal heterostructures 
[39]. Importantly, in such multi-orbital models, topological 
order can emerge from purely local interactions, whereas 
nonlocal interactions are typically required for single- 
orbital models [39]. Another possible setting to observe 
interaction-generated topological states is stacked graphene 



[40,41]. Finally, three-dimensional topological insulators 
can also be generated from interactions [133]. 

Topological states arising from electronic interactions 
have also been discussed in the context of systems with 
intrinsic spin-orbit coupling that undergo a Mott transition 
from a regular to an exotic quantum spin Hall phase of 
spinons, in which interaction-generated low-energy spin 
excitations inherit the topological band structure [29,63, 
108, 134], see also section 4.3. Alternatively, interactions 
can give rise to a transition from a nontopological state 
to an integer quantum spin Hall state, as in [73,84], 
essentially via a renormalization of the effective spin-orbit 
coupling. In relation to fractional states to be discussed 
next, it is important to distinguish fractionalized and 
fractional topological insulators. The former are states 
with spin-charge separation but typically integer charge 
excitations, whereas the latter are states that exhibit charge 
fractionalization as in the fractional quantum Hall effect. 

4.5. Fractional topological insulators 

Given the close relation between integer quantum Hall 
states and quantum spin Hall insulators, Bernevig and 
Zhang [135] envisaged the possibility of fractional time- 
reversal invariant insulators arising from strong electronic 
interactions [135]. Levin and Stern [32] considered the case 
of independent fractional quantum Hall states for j" and l 
spins, and established that the criterion for the existence 
of edge states protected by time reversal symmetry is 
that cr®/e* is odd, where e* is the elementary charge 
of the fractional Hall state. However, a combination of 
two fractional quantum Hall states is not necessarily a 
topological insulator with protected edge modes, but can 
instead also lead to fractional trivial insulators with o s xy /e* 
even [32]. In the latter, edge modes can be gapped out by 
time-reversal invariant perturbations. 

A lattice model that supports the integer quantum Hall 
effect without a magnetic field was introduced by Haldane 
[8], see (5). Analogous models for fractional quantum 
Hall states were discovered only recently, for example, in 
[30,81, 136-140]. A common feature of such models is 
the existence of flat bands with a nonzero Chern number 
(so-called Chern bands) that mimic the Landau levels of 
more traditional fractional quantum Hall states observed in 
strong magnetic fields [2]. The flatness allows electronic 
interactions to be large compared to the kinetic energy 
but small compared to the band gap (thereby avoiding 
level mixing effects and allowing for a projection onto a 
single-band problem). A number of Haldane-like models 
with these properties where suggested [30, 123, 136-138], 
and the existence of a fractional quantum Hall phase was 
demonstrated in [30, 138]. The prospect of fractional 
quantum Hall states without magnetic fields has also 
inspired the search for materials with partially filled Chern 
bands arising from spontaneous orbital order, see [31] 
and references therein. Similar to the case of topological 
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Mott insulators (section 4.4), the necessary spin-orbit 
coupling (or spin-dependent magnetic flux) can arise from 
interaction-induced spontaneous symmetry breaking [31]. 
An adiabatic connection between Chern insulators and the 
Hofstadter problem is discussed in [141]. 

In the integer quantum Hall effect, the chemical 
potential lies in the gap between the highest filled and 
the lowest unfilled energy levels. As long as the energy 
scale for electronic interactions is smaller than the gap 
to excited states, interaction effects can be considered as 
a small perturbation. In contrast, electronic correlations 
are at the heart of fractional quantum Hall liquids and 
hence also fractional topological insulators; for a review 
see [2]. At the special filling fractions for which a xy 
exhibits a plateau, a family of ground states with well- 
defined topological degeneracy is spontaneously selected 
from the highly degenerate set of many-body states. This 
process is inherently nonperturbative, and the fractional 
state has quasiparticle excitations with fractional charge 
and fractional statistics. The extremely rich physics 
known from the fractional quantum Hall effect stimulates 
research on fractional topological insulators. A previously 
encountered example of such a fractional spin Hall liquid is 
the QSH* phase of the SI model [76], see section 4.3. 

Neupert et al [98] considered the general case of a 
fractional liquid with time reversal symmetry but no spin 
conservation. Such states (including the noninteracting, 
integer quantum spin Hall states) are described at low 
energies by Abelian Chern-Simons theories, whereas the 
special case with spin conservation is described by a 
U(l) x U(l) Chern-Simons or BF theory [53, 142]. If 
spin is not conserved, the existence of protected edge 
states can be related to a Z 2 invariant independent of 
at. that reduces to the usual Z„ index in the integer, 

xy z 

noninteracting limit [98]. A general Chern-Simons theory 
of two-dimensional, time-reversal invariant insulators with 
Abelian quasiparticles and a classification of fractional 
topological insulators was given by Levin and Stern [143]. 
Recent work along these lines includes parton [144] and 
composite fermion [145] descriptions, the relation between 
fractional topological insulators and quantum Hall states 
in Landau levels [146], and an effective field theory of 
general fractional two-dimensional topological insulators 
with spin-orbit coupling [147, 148]. Because the fractional 
states arise from strong interactions, topological order will 
typically compete with other many-body instabilities such 
as charge-density-wave formation. For Chern insulators, 
the stability of the topological state has been found to 
be enhanced by a certain amount of band dispersion or 
for bands with higher Chern numbers [149]. Finally, the 
experimental detection of fractional topological insulators, 
including the determination of the filling fraction from 
transport measurements, is discussed in [150]. 

An important open problem is the question under 
what circumstances two coupled fractional Hall liquids 



form a fractional topological insulator with protected edge 
modes. The interactions between the two fractional states 
may completely destroy the insulating state, give rise 
to spontaneous breaking of time-reversal symmetry, or 
lead to novel, non- Abelian states [98]. It will also be 
interesting to see if the physics beyond the Laughlin states 
recently reported for Chern insulators [140, 151] carries 
over to fractional topological insulators. Most of the 
numerical work on fractional topological insulators so far 
relies on exact diagonalization, with an exponential scaling 
of computational effort with system size. A possible 
alternative is the use of tensor network states [152]. 

5. Edge correlation effects 

The existence of metallic, protected edge states is one of the 
most obvious manifestations of the topological properties of 
quantum spin Hall systems [9, 10], and has played a crucial 
role for the experimental verification of the existence of 
topological insulators [14, 15]. The fundamental properties 
of these edge states can be understood in the noninteracting 
limit, as reviewed in section 3. Due to time-reversal 
symmetry, the edge states always exist as Kramers pairs, 
and electrons of opposite spin propagate in opposite 
directions along a given sample edge. For an odd number 
of such helical edge state pairs, time-reversal symmetry 
provides a protected, gapless crossing point [10]. In 
contrast, single-particle backscattering is possible for an 
even number of pairs, and the edge states are not protected 
[9, 10]. The parity of the number of edge state pairs can 
therefore be related to the Z 2 topological invariant [9]. 

Compared to the gapped bulk, electron-electron 
interactions can be expected to strongly affect the gapless 
edge states [55]. Most importantly, in the presence of 
interactions, inelastic scattering becomes possible, and 
is not forbidden by Kramers theorem, as it connects 
states of different energy [55]. Moreover, whereas time- 
reversal symmetry forbids elastic single-particle processes, 
two-particle backscattering is allowed even for a single 
Kramers pair. The possibility of such processes, which 
can spontaneously break time-reversal symmetry at the 
edge [54,55], is a crucial difference to the quantum 
Hall effect, where counter-propagating states (and hence 
backscattering) are completely absent [94]. As pointed 
out in [54,55], the qualitative, topological distinction 
between odd and even numbers of edge state pairs in 
the noninteracting case becomes a quantitative distinction 
(depending on the interaction strength) in the presence 
of correlations. However, as discussed below, for 
the class of models with conserved spin and without 
impurities, the topological protection and the bulk- 
boundary correspondence remain intact. 

An odd number of Kramers edge state pairs, and 
hence a topologically protected quantum spin Hall state, 
can be realized only at the boundary of a two-dimensional 
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system [54], giving rise to the notion of a holographic 
liquid [54]. Hence, edge states of a quantum spin Hall 
insulator may be regarded and theoretically described (to 
a good approximation) as one-dimensional liquids that 
inherit their topological properties from the bulk [55,95]. 
Similar to other quasi-one-dimensional systems such as 
quantum wires, quantum fluctuations play a key role in 
determining the physics of edge states. Consequently, 
whereas mean-field theory provided valuable insight into 
the effect of bulk correlation effects, it is insufficient for 
edge physics and in fact gives misleading results as a 
result of instabilities [111]. Similarly, given the non-Fermi 
liquid character of the metallic edge states, perturbative 
approaches fail. However, Luttinger liquid theory [153] 
provides a rather accurate description of many aspects [54, 
55]. Although a strictly one-dimensional description seems 
justified by the exponential suppression of scattering to 
bulk states, signatures of the holographic nature of the 
helical edge states are visible as quantitative deviations 
from bosonization predictions [113]. 

In this section, an introduction to the helical Luttinger 
model is given [54,55] (for a previous review see [12]), 
followed by a discussion of analytical and numerical results 
with an emphasis on theoretical aspects of the intensely 
studied integer quantum spin Hall insulators. A more 
experimentally oriented review of transport properties was 
given by Tkachov and Hankiewicz [154]. Other interesting 
developments not addressed in the following include 
Majorana edge states of topological superconductors [116], 
and edge states of fractional topological insulators [32, 
150]. Finally, the stability of edge states in magnetically 
ordered topological insulators is discussed in [56, 1 14, 121— 
123]. 

5.1. Luttinger liquid description 

The low-energy description of helical edge states in terms 
of Luttinger liquid theory was established in two early 
papers by Wu et al [54] and Xu and Moore [55]. The 
corresponding theory and hence the physics is independent 
of the particular model for the bulk as long there is an 
adiabatic connection to the noninteracting case. (For edge 
state theories of fractional topological insulators see [32, 
98, 150].) However, details such as the number of edge 
states, the allowed interactions and the location of the 
crossings in momentum space are determined by the Z 2 
topological invariant, the spin symmetry, the bulk band 
structure and the edge topology (for example, zigzag or 
armchair for honeycomb lattice models). For a generic 
topological insulator, the low-energy model for the edge 
states can be derived from Chern-Simons theory [55,95], 
see also section 3. 

In the following, a single pair of edge states is 
considered, assuming a geometry that is wide enough for 
the two edges to be independent. Correlation effects 
on even numbers of edge state pairs, corresponding to a 



topologically trivial state in the noninteracting case, have 
been discussed by Xu and Moore [55]. Surprisingly, in 
this situation, the edge states can be stabilized against 
localization by interactions [55]. Effects of inter-edge 
tunnelling [74, 111] will be addressed in section 5.2. 

The following discussion follows Wu et al [54] (see 
also [12]). Assuming a system with periodic boundaries in 
the x direction and open boundaries in the y direction, and 
linearizing the edge dispersion as e (T (fc) = ±v F (k — k F ), 
the free continuum theory takes the form 

H Q =v F J Ax (VRt^Rt - ^14^14) • (30) 

The helicity of the edge states is apparent from the existence 
of right moving electrons (R) with spin j" and velocity v F , 
and left moving electrons (L) with spin J, and velocity 
—v F . The spin directions can be defined with respect to 
an arbitrary quantization axis even if spin is not conserved. 
The Fermi velocity v F is related to the bulk band gap and 
the edge topology, but analytical expressions for specific 
models are in general not known (see, for example, [155]). 
Because of the correspondence between the direction of 
motion and spin, a spin index is not required, and drops 
out of the effective low-energy Luttinger liquid theory (see 
below). An interesting consequence of the reduced number 
of degrees of freedom compared to ordinary Luttinger 
liquids is the predicted possibility to observe fractional 
charges [156]. 

In constructing the low-energy model, it is important 
to consider the possible interactions. Elastic perturbations 
which are odd under time reversal, most importantly 
single-particle backscattering, are not allowed. Inelastic 
single-particle scattering is allowed, but requires a spin 
flip. It can arise, for example, from magnetic impurities, 
Rashha coupling in combination with phonons [157], or a 
coupling to spin-1 bosons (for example, spin fluctuations 
[113]). Forward scattering as well as elastic or inelastic 
two-particle processes are generally allowed [54,55]. 
Although a discussion of scattering processes based on 
the single-particle picture is physically intuitive, a many- 
body formulation [55] is more appropriate for interacting 
systems, and leads to the same low-energy Luttinger 
model. A continuum theory of helical edge states was 
presented in [158]. Bosonization results for non-helical 
one-dimensional systems with spin-orbit coupling can be 
found in [159, 160]. 

Forward scattering is described by [54] 

H fw = 9 J dx ^R^Rt^U^U ■ (31) 

In clean systems, the most important two-particle term is 
umklapp scattering (i.e., backscattering from a periodic 
lattice potential), given by [54] 

H um = 9u J dxe-^^(x)^(x + a) 

x^ u (z + a)V> u (a:)+H.c. (32) 
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Whereas H iw is generically present in helical liquids, the 
umklapp term is rather special as it requires half filling 
(4/c F = 2wl, I integer) and is only generated in the 
presence of spin-nonconserving terms in the Hamiltonian, 
such as the Rashba spin-orbit interaction [161]. It is 
therefore generically absent in models with U(l) spin 
symmetry, even at half filling. Two-particle backscattering 
can also arise from magnetic impurities [54, 162-164] or 
quenched disorder [10,54], see also section 5.2. Finally, 
in systems without spin conservation, combined single- 
particle forward- and backscattering can become important 
[71]. 

Under the assumption that the edge states remain 
quasi-one-dimensional, interaction effects can be studied 
in the powerful framework of bosonization. If only 
forward scattering is allowed, electron-electron interactions 
can be accounted for by the Luttinger liquid parameters 
K = y/(v F - g)/(v F + g) and v = ^uf - g 2 (these 
expressions are valid for g <C v F ). The corresponding 
helical Luttinger model reads 



1.00 



K{djf 



(33) 



where <ft = R + (f> L and 9 = R — </> L ; the possibility of 
using these variables to obtain a spinless model is a direct 
consequence of helicity. 

Because umklapp- and backscattering are excluded, 
the quadratic and hence exactly solvable Luttinger 
model (33) describes correlated but metallic helical edge 
states, as realized in a number of spin-conserving theoreti- 
cal models without disorder. Electron-electron interactions 
lead to a change of the parameters K and v: K = 1 and 
v = v F for noninteracting electrons, and K < 1 (K > 1) 
for repulsive (attractive) interactions. However, no quantum 
phase transitions occur, and all correlation functions decay 
with power laws with exponents determined by K. Explic- 
itly, the equal-time correlation functions can be derived as 

cos(2fc F r) , 



S xx {r) = (S x {r)S x (0)} 
S zz {r) = (S z (r)S z (0)) 
= (n(r)n(0)) ~ 



(34) 



1 



N(r) 
P(r) 



Rt 



(r)^(r)4. t (0)^(0)) 



r2/K 



corresponding to transverse spin, longitudinal spin, charge 
and pairing correlations, respectively. The results 
in (34) imply that transverse spin correlations (described 
by S xx (r)) are the slowest decaying correlations for 
repulsive interactions, whereas pairing correlations P(r) 
dominate for attractive interactions. The longitudinal 
spin correlations S xx (r) (here S zz ^ S xx reflects the 
broken SU(2) spin symmetry in the presence of spin-orbit 
coupling) and charge correlations N(r) decay exactly as in 
a Fermi liquid. Finally, only S xx (r) depends on the Fermi 
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Figure 13. (a) Transverse spin correlator for different interaction strengths 
U/X (A: spin-orbit coupling), showing a power-law decay with an 
interaction-dependent exponent in accordance with (34). Lines are fits 
to the form A/r ,rl . Data are normalized to 1 at r = 1 for comparison. 

(b) Renormalized transverse spin structure factor at q = 0, showing the 
absence of long-range magnetic order in the thermodynamic limit. Lines 
are fits to the form S xx (q = 0)/L = b/ L + c/ L 7 ) , with r] taken from (a). 

(c) Transverse spin structure factor at U/X = 8 for different band fillings 
n. All results were obtained from zero-temperature quantum Monte Carlo 
simulations of an effective model with Hubbard interactions only at the 
edge of a zigzag ribbon. Data taken from [113]. 



momentum via the term cos(2fc p r), which can be related to 
spin-flip processes [113]. 

The bosonization predictions (34) have been tested 
numerically by considering the quantum spin Hall phase 
of microscopic, interacting models. To verify the 
characteristic interaction-dependent power-law decay of 
transverse spin and pairing correlation functions, a hallmark 
feature of one-dimensional liquids predicted by (34), rather 
large system sizes are required. Existing work along these 
lines is based on quantum Monte Carlo simulations of 
the KMH model [103], and of an effective model that 
takes into account electron correlations only at the edge 
[99, 113]. This approach is justified by the minor role of 
bulk interactions in the topological band insulator phase, as 
discussed in section 4. The effective model allows to access 
large system sizes without any further approximations. 
(The results shown in [99,113] were obtained with an 
incorrectly implemented bulk Green function, and therefore 
correspond to an alternative model of a quantum spin Hall 
insulator rather than to the KM model. However, using 
the KM Green function leads to very similar results, and 
none of the conclusions are affected, see errata of [99, 



CONTENTS 



21 



113]. In the following, figure 15 shows results for the KM 
model with edge interactions, whereas figures 13, 14, 18 
correspond to previously obtained results for the alternative 
model.) 

Results for S xx (r) are reproduced in figure 13(a) 
[113]. For a local Hubbard interaction, the forward 
scattering matrix element g depends on U, whereas the 
Fermi velocity is determined by the spin-orbit gap and 
hence by A go . The effective interaction strength for 
the edge states is therefore set by the ratio U/X so 
[113]. Figure 13(a) reveals that transverse spin correlations 
S xx (r) decay increasingly slower with increasing U/A so , 
corresponding to a reduction of K in (34). All other 
correlation functions decay much faster [113]. Similar 
results have been obtained in a full bulk calculation of the 
KMH model [103, 104]. Despite the strongly enhanced 
transverse spin correlations, the Luttinger model predicts 
the absence of long-range order. This prediction is 
consistent with the finite-size scaling of the spin structure 
factor shown in figure 13(b). Figure 13(c) confirms the 
doping dependence of S xx (r) by showing its Fourier 
transform (the spin structure factor) at different particle 
densities. Away from half filling (n — 1, k F = it), 
the cos(2fc p r) factor in (34) becomes visible. From the 
numerical results, the Luttinger parameters K and v can be 
determined using finite-size scaling. The results in figure 14 
demonstrate a strong suppression of K with respect to the 
noninteracting value K = 1, and a slight renormalization of 
the velocity in agreement with [105, 106]. The parameter 
K has also been determined from quantum Monte Carlo 
simulations of the KMH model by fitting power laws to 
correlation functions at a fixed system size [103], see 
figure 17. The renormalization of the velocity v was also 
studied in [106]. 

While the numerical value of K for given microscopic 
parameters U and A go will depend on the details of the 
model, repulsive interactions generically lead to dominant 
2fc F transverse spin fluctuations on helical edges [54, 55, 
99, 103, 111,113]. In the light of the bulk results discussed 
in section 4, transverse spin correlations are therefore a 
hallmark feature of correlated topological insulators. This 
finding is expected to hold even in the presence of Rashba 
coupling, as long as interactions are not too strong (see 
below). Although mean-field theory is not appropriate to 
describe correlated edge states, it does provide some insight 
into the differences between bulk and edge correlations. 
For a system with periodic boundaries, mean-field theory 
gives a finite critical value U c at which the topological 
band insulator undergoes a quantum phase transition to a 
magnetic Mott insulator [63], see inset of figure 10. In 
contrast, a system with edges is unstable toward the opening 
of a gap due to a nesting-related, logarithmically diverging 
susceptibility, and shows a gapped magnetic state with 
long-range order, (S xx ) ^ 0, for any finite U [111]. This 
so-called edge magnetism has been studied in the context 
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Figure 14. Luttinger liquid parameters K and v for an effective model 
of helical edge states, as obtained from quantum Monte Carlo simulations 
and finite-size extrapolation. Data taken from [113]. 

of graphene [165-168]. Lee [111] has shown that for weak 
enough interactions, taking into account fluctuations around 
the mean-field saddle point restores the expected power-law 
decay at the edges. Although long-range order is absent, 
exact numerical results for the KMH model expose that 
magnetic correlations are still strongest at the edge, and 
decay quickly with the distance from the edge [103]. 

The effective model with interactions only at the edge 
sites provides an accurate description of the edge states 
inside the quantum spin Hall phase [99, 113]. In addition 
to capturing the power-law decay of correlations expected 
from Luttinger liquid theory, this approach also permits to 
calculate dynamical correlation functions on large systems. 
Figure 15(a) shows the single-particle spectral function 
of the KMH model with edge Hubbard interactions on a 
zigzag ribbon at A go /i = 0.2 and U/t = 2. As in the 
noninteracting case (figure 4(a)), there is a gapless crossing 
point at k = ir. The absence of any features related to 
spin-charge separation is due to helicity, as reflected in the 
spinless effective model (33), and is a generic feature of 
one-dimensional systems with spin-orbit coupling [159]. 
The single-particle spectrum of a correlated helical liquid 
has also been obtained with cluster dynamical mean-field 
theory on zigzag and armchair edges [105], see figure 16, 
for the KMH model using the variational cluster approach 
[106], and for the BHZH model [71,73,74]. The real- 
space Green function of the helical Luttinger model has 
been calculated analytically [169]. 

Whereas dynamical two-particle correlation functions 
are not reliably accessible with cluster methods [105, 106] 
and difficult to obtain within bosonization, they can be 
calculated with high momentum resolution using quantum 
Monte Carlo simulations of the effective model [99, 113]. 
Figure 15 shows results for the KMH model with edge 
Hubbard interactions. The charge structure factor in 
figure 15(b) and the longitudinal spin structure factor in 
figure 15(c) are characterized by a linear mode at long 
wavelengths that has the same velocity as the edge state 
visible in the single-particle spectrum shown in figure 15(a). 
The corresponding spectral weight is much smaller in 
the charge channel than in the spin channel. The linear 
mode is related to spin-conserving excitations within the 
helical states, whereas excitations at higher energies involve 
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bulk states. The transverse spin structure factor shown in 
figure 15(d) involves spin flips — connecting the two helical 
edge modes — and therefore exhibits a continuum of low- 
energy excitations near q = [113]. 

As confirmed by the explicit numerical calculations 
discussed in this section, for weak enough interactions, the 
fixed point described by the forward scattering model (33) 
qualitatively describes the physics of helical edges. As long 
as K does not fall below the respective critical values to be 
discussed below, the model remains valid even for systems 
with nonconserved spin [71], disorder [10, 54], or magnetic 
impurities [54, 162]. 

5.2. Strongly correlated regime 

The robustness of helical edge states with respect to 
interactions is a result of time-reversal invariance. At 
the time-reversal invariant momenta k = 0, it, the two 
edge modes necessarily have the same energy, e^(fc) = 
e^(fc), so that the system is gapless. A gap in the 
edge states can arise from a breaking of time-reversal 
invariance either globally (for example, due to spontaneous 
symmetry breaking and the onset of magnetic order in the 
bulk, see section 4), or only at the edge. In the latter 
case, the bulk remains a topological insulator, so that the 
bulk-boundary correspondence breaks down [93, 170, 171]. 
This section explores the conditions for the occurrence of 
gapped edge states, and reviews the consequences of strong 
interactions — compared to the bulk gap — in systems that 
remain gapless for symmetry reasons. In general, a gap 
in the edge states requires strong interactions, as measured 
by the Luttinger liquid parameter K. Evidence for strong 
interactions (K <C 1) for a wide range of microscopic 
model parameters comes from numerical simulations [103, 
113], see figures 14 and 17. 

The impact of strong interactions can be analyzed 
by combining bosonization with the renormalization group 
method [54, 55]. Whereas forward scattering does not open 
a gap, backscattering terms become relevant for strong 
enough interactions [54,55]. The umklapp term, important 
at half filling and for systems without spin conservation, 
becomes a relevant perturbation for K < 1/2 [54,55, 
111]. Backscattering due to quenched disorder leads 
to a gap if K < 3/8 [54,55], whereas backscattering 
involving magnetic impurities becomes relevant for K < 
1/4 [162,163]. Finally, a gap in the edge states can 
also open as a result of inter-edge scattering, relevant for 
AttK + (AttK)^ 1 < 4 according to Lee [111], and in the 
presence of a lattice of magnetic impurities when K < 
1/2 [163]. The helical liquid hence shows remarkable 
differences compared to ordinary Luttinger liquids. In 
the latter, if K < 1, transport is completely blocked at 
T = by either disorder [96] or a single magnetic impurity 
(the Kondo problem) [172]. In contrast, a helical liquid 
is unconditionally stable with respect to single-particle 
backscattering [10], and can avoid the site with the Kondo 




Figure 15. Dynamical correlation functions for (a) single-particle, (b) 
charge, (c) longitudinal spin, and (d) transverse spin excitations. Results 
obtained from zero-temperature quantum Monte Carlo simulations of the 
KMH model on a 24 X 64 zigzag ribbon with a Hubbard U only at the 
edge sites [99, 113]. Here U/t = 2 and X so /t = 0.2. 



singlet by penetrating into the bulk provided K > 1/4 
[162, 163]. The observation of metallic edge states in 
experiments on HgTe quantum wells [14-16], in which 
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disorder and impurities may be expected to be present to 
some degree, suggests that interactions are rather weak, so 
that the simple model (33) is sufficient. 

Because of the complexity of the many-body problems 
with either additional Rashba interaction, impurities or 
disorder, the above predictions have not yet been verified 
by simulations of microscopic models. Nevertheless, 
significant progress in the understanding of correlated 
helical edges has been made. The simplest scenario 
involves simultaneous breaking of time-reversal symmetry 
in the bulk and at the edge, as realized, for example, 
in the KMH model at large U/t (section 4.3). Because 
time-reversal symmetry is broken in the magnetic phase, 
a topologically trivial state without gapless edge states is 
expected. The evolution of the edge states across the 
magnetic bulk transition can only be studied in a model 
with an interacting bulk. Hence, the effective model of 
[99, 113] is not applicable. Instead, this problem has been 
studied in the framework of the KMH model by Wu et 
al [105] using cluster dynamical mean-field theory, and 
by Yu et al [106] using the variational cluster approach. 
Although these results do not fully capture the Luttinger 
liquid properties of the edge states, the results agree well 
with expectations. Figure 16 shows the single-particle 
spectral function for an armchair edge as obtained by Wu 
et al [105]. For the value A so /t = 0.2 considered, the 
critical value for the bulk magnetic transition within their 
approximation is U/t ~ 5.2, which seems to agree with 
the opening of a gap in the edge states [105]. The low- 
energy model for the gapped edge states can be written 
as -ff HL — 9 m sin(v / 47T0) (cf. (33)), where m is the 
bulk magnetization that, if nonzero, leads to a relevant 
perturbation for repulsive interactions [105]. A mean-field 
analysis with additional spin-wave fluctuations has been 
given by Lee [1 11]. 

To observe a topological insulator with Mott- 
insulating edges, time-reversal symmetry has to be broken 
at the edge. The observed tendency toward strong magnetic 
correlations follows already from the simple model (33), 
which predicts long-ranged but power-law transverse spin 
correlations for K -C 1, see (34). The exponent K has 
been calculated numerically for different (spin-conserving) 
models [103, 113], and it was found that it can be substan- 
tially smaller than 1, see figures 14 and 17. However, as 
shown in figure 13, and in accordance with Luttinger liquid 
theory, true long-range order is absent even for very strong 
interactions [113]. Assuming that the edge states can be 
considered to be strictly one-dimensional, the absence of 
long-range order is implied by the Mermin- Wagner theo- 
rem, which forbids continuous symmetry breaking in one 
dimension. However, the continuous U(l) spin symme- 
try can be reduced to a discrete Z 2 Ising symmetry in the 
presence of Rashba spin-orbit coupling, or in more general 
models such as (9). In this case, umklapp processes of the 
form (32) become possible, and Ising long-range magnetic 
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Figure 16. Single-particle spectral function of the KMH model on an 
armchair ribbon [105] with spin-orbit coupling A = Ago an d Hubbard 
interaction U . Results were obtained with the cluster dynamical mean- 
field theory. The top row corresponds to the quantum spin Hall phase, 
whereas the bottom row shows the spectrum in the magnetic (left) and the 
nonmagnetic insulator (right), respectively. (Reprinted with permission 
from [105]. Copyright 2012 by the American Physical Society). 

order becomes allowed at T = if K < 1/2 [54,55, 111]. 
However, at T > 0, time-reversal symmetry will be restored 
by thermal fluctuations. This mechanism permits the exis- 
tence of a topologically nontrivial quantum spin Hall in- 
sulator with gapped edge states and time-reversal symme- 
try, and causes the bulk-boundary correspondence to break 
down [93, 170, 171]. Temperature-dependent edge and bulk 
magnetism in a quantum spin Hall insulator has been dis- 
cussed by Shitade et al [22]. The possibility of umklapp- 
driven magnetic order at the edge has been taken into ac- 
count by Zheng et al [103] in their phase diagram of the 
KMH model. The latter, shown in figure 17, depicts a rather 
large region where K < 1 /2, although the simulations were 
done for A R = 0. (The gap in the edge Green functions in 
[103] is a finite-size effect.) 

Based on the above arguments, and for models with 
conserved spin, the edge states are expected to become 
gapped exactly at the point where long-range magnetic 
order in the bulk breaks time-reversal invariance and hence 
destroys the topological state. A slight complication is 
that the helical liquid is holographic, and hence exists 
only at the boundary of a two-dimensional system [54], 
whereas the absence of continuous symmetry breaking 
applies only to strictly one-dimensional systems. The two- 
dimensional bulk is not taken into account in the Luttinger 
model (33), but is included in the effective model with 
interactions only at the edge; no evidence for long-range 
order has been found (see figure 13) [113]. For graphene 
edges, it has been argued that the Mermin- Wagner theorem 
is invalidated by the momentum dependence of electron- 
electron interactions, and numerical results are consistent 
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Figure 17. Phase diagram of the half-filled KMH model with spin-orbit 
coupling A = Ago and Hubbard interaction U from quantum Monte Carlo 
simulations [103]. Also shown are regions where the Luttinger parameter 
K < 1/2 and K > 1/2, respectively. (Reprinted with permission from 
[103]. Copyright 201 1 by the American Physical Society). 

with long-range magnetic order [167, 168]. 

The difficulties associated with capturing large length 
scales (cf. figure 13) and quantum fluctuations, as well 
as with the unambiguous detection of topological states, 
have led to several claims for the existence of edge 
Mott insulators (or edge superconductors) in systems with 
conserved spin. Paramagnetic, inhomogeneous slave-rotor 
calculations for the KMH and the BHZH model suggest the 
existence of a paramagnetic edge-Mott state for large values 
of A go [70]. Similar conclusions have been drawn based 
on variational Monte Carlo simulations of the KMH model 
[104]. Finally, a symmetry -breaking superconducting state 
has been reported for the KMH model with attractive 
interactions [173] using a mean-field calculation. In two 
dimensions, the slave-rotor result [70] is known to be 
unstable with respect to gauge fluctuations [63]. A similar 
but stable phase (a topological Mott insulator) has been 
predicted based on slave-boson calculations [108]. The 
variational Monte Carlo results for the spin and charge 
Drude weights [104] have rather large error bars, and the 
authors interpret them either as a transition or a crossover. 
The latter possibility agrees with other numerical work [99, 
105, 113], see discussion below. 

The existence of a superconducting edge state can be 
excluded by exploiting the exact particle-hole symmetry of 
the KMH Hamiltonian, and the canonical transformation 
c^j. H> (— l) l c,.^, with (— l) 1 = ±1 depending on 
the sublattice. Similar to the Hubbard model, this 
transformation exchanges the role of spin- 2 and charge 
operators in the Hamiltonian, and flips the sign of the 
Hubbard term [174]. The dominant transverse spin 
correlations observed in the repulsive model translate into 
dominant pairing correlations P(r) in the attractive model. 
However, all correlation functions decay with a power law, 
as can be explicitly confirmed by quantum Monte Carlo 
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Figure 18. Luttinger liquid parameter K and spectral weight of long- 
wavelength charge excitations Z from quantum Monte Carlo simulations 
of an effective model for helical edge states. (Reprinted with permission 
from [113]. Copyright 2012 by the American Physical Society). 

simulations for the attractive case. The symmetry breaking 
observed in [173] at the mean-field level therefore has the 
same origin as the long-range magnetic order in repulsive 
model, namely a logarithmically divergent susceptibility 
[111]. 

Although no gap is opened in quantum spin Hall 
insulators with conserved spin, strong interactions lead 
to important modifications of the edge states. Apart 
from a renormalization of the velocity v with increasing 
l7/A so that follows from the low-energy model (33), see 
figure 14, several numerical studies indicate a noticeable 
suppression of charge transport, visible for example in the 
dynamical charge structure factor at long wavelengths and 
low energies in figure 15(b), or in the Drude weight [104], 
with a simultaneous enhancement of spectral weight for 
transverse spin excitations [99,104,113]. Similarly, the 
spectral weight of the edge states in the single-particle 
spectrum is diminished by strong interactions [99, 105, 
113]. 

The Drude weight, D, can be related to the spectral 
weight, Z, of the linear mode in the dynamical charge 
structure factor via D = Zv (v being the velocity); Z 
can be measured in quantum Monte Carlo simulations 
[113], and results are shown in figure 18. On the other 
hand, for a Luttinger liquid described by (33), the Drude 
weight is given by D = Kv [113], so that the apparent 
difference between Z and K in figure 18 indicates that 
the model (33) is incomplete. While Z = K in the 
noninteracting limit, the Drude weight is significantly 
smaller than predicted theoretically (i.e., smaller than D = 
Kv) for nonzero interactions. The discrepancy has been 
attributed to inelastic spin-flip scattering mediated by the 
pronounced magnetic fluctuations that exist for K -C 1, 
and to scattering processes that involve bulk states [113]. 
Whereas the importance of bulk states can be expected if 
the interaction U becomes comparable or larger than the 
bulk band gap, the numerical results suggest that deviations 
occur even for rather weak interactions [113]. 
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6. Topological invariants of correlated systems 

The discovery of the remarkable relation between the 
quantized Hall conductivity and the first Chern number 
in the context of the quantum Hall effect [3-5] has 
played a key role for the understanding of topological 
phases of matter. In section 3, it has been shown that a 
topological invariant can also be defined for quantum spin 
Hall insulators, and that the latter is — in the simplest case — 
determined by the difference of the Chern numbers of 
two quantum Hall states with opposite Hall conductivities 
[62]. Reviews of topological invariants with a focus on 
noninteracting systems can be found in [12,20,36]. Here 
the emphasis is on the calculation of topological invariants 
for correlated systems by means of numerical methods. 
More general discussions can be found, for example, in 
[20, 175]. 

The direct calculation of the topological invariant of 
interacting systems is in general rather difficult. Instead, 
the topological properties can in certain cases be inferred 
indirectly either via the bulk-boundary correspondence, 
or based on an adiabatic connection to a noninteracting 
system. 

Because time-reversal symmetry protects topological 
insulators against weak interactions, the bulk-boundary 
correspondence can in principle be used to determine 
whether or not a topological state survives upon switching 
on electron-electron interactions. However, this approach, 
followed for example in [73, 105, 106], faces a number 
of difficulties. In addition to the inherent difficulty of 
obtaining numerical results on large enough systems with 
open boundaries to demonstrate the existence of gapless 
excitations, correlation effects and electron tunnelling 
between edges can gap out the edge states despite the 
existence of topological order in the bulk (see section 5.2). 
Although gapped edge states can be reconciled with a 
nonzero topological invariant theoretically [93, 170, 171], a 
distinction of topological and trivial states based on edge 
states alone will be difficult in actual numerical calculations 
for strongly correlated systems. Nevertheless, the bulk- 
boundary correspondence is a valuable theoretical concept. 
For example, it has been used to formulate criteria for the 
stability [32,98] and experimental identification [150] of 
fractional topological insulators. 

A topological phase may also be identified by 
exploiting an adiabatic connection to a noninteracting state 
(for which the invariant can easily be calculated). For 
example, as discussed in section 4.2, the KMH model 
with A R = describes a quantum spin Hall state for 
U = and A so ^ 0. For small to intermediate values 
of U ft, the system remains fully gapped and time-reversal 
invariant, and can be adiabatically tuned (via U/t) into 
the noninteracting quantum spin Hall state without closing 
the band gap. As emphasized in section 4.5, no such 
adiabatic connection to a noninteracting state exists, for 
example, in the case of fractional topological insulators or 



topological Mott insulators. Finally, even with the use of 
exact numerical methods, it is not always possible to decide 
if an excitation gap closes or not across a transition [100]. 

Given the limitations of these indirect probes for 
topological states, it is important to be able to directly 
calculate the topological invariant. For spin-conserved 
quantum spin Hall insulators, it is sufficient to determine 
the Chern number for one spin sector. Because the Chern 
number is related to the quantized Hall conductivity, it can 
in principle be obtained without approximations from the 
Kubo formula [87, 88]. The Z 2 invariant then follows from 
(15). Nevertheless, for correlated electron systems, even 
the calculation of the conductivity can be quite demanding. 
When spin is not conserved, the separation into two spin 
sectors is obviously no longer possible, and other concepts 
such as the spin Chern number are required. In particular, 
even an exact calculation of the Hall conductivity via the 
Kubo formula is insufficient because a* is no longer 
quantized [62]. 

Considering first the integer, spin conserved case, (13) 
and (14) are not applicable to interacting systems, because 
they rely on a representation in terms of Bloch states. 
Instead, the Chern number in the presence of interactions 
and/or disorder can be calculated from [87] 
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Compared to (13), the Bloch states \u m ) are replaced by 
the many-body ground state |\& ), and the momenta k , 

(j> describing 



k y are replaced by phase parameters . 
general boundary conditions \I> (. . . , {x 
gKPxijr f j {a; y.}, . . .) (and similar for the y direction) 
[87]. On a torus geometry, the phase shifts <f> x , (f> y can be 
related to magnetic fluxes [88]. Assuming the existence 
of a bulk energy gap and a unique, nondegenerate ground 
state, C as defined by (35) is an integer Chern number 
directly related to the Hall conductivity [87]. An alternative 
representation is given by [87, 88] 
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The derivatives of the Hamiltonian are numerically easier 
to evaluate than derivatives of the wave functions [88], 
and (36) reveals the relation to the Kubo formula [88]. For 
simplicity, the zero-temperature limit has been taken. At 
finite temperatures, the Hall and spin Hall conductivities are 
in general not exactly quantized, see for example [72, 176]. 

Equations (35) and (36) can be extended to the case 
of a fractional quantum Hall state with a d-dimensional 
ground-state manifold {|*^)}, K = l,...,d. The Chern 
number reads [87] 
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The corresponding generalization of (36) to the fractional 
case can be found, for example, in [177]. Equations (35)- 
(37) are particularly useful in combination with exact 
diagonalization methods, and recent applications in the 
context of correlated Chern insulators include [31,66,68, 
177]. Twisted boundary conditions can also be used in 
combination with quantum Monte Carlo simulations [66]. 

A generalization for integer quantum spin Hall 
insulators without spin conservation is based on the Chern 
number matrix [62, 91] 
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with spin-dependent boundary twists 4> x , <j) y ■ If s P m i s 
conserved, G^ = = whereas = = 

±1. Sheng et al [62] defined a total charge Chern number 
C c = X^ercr' C' 7 ' 7 an d a tota l s P m Chern number C s = 
^2 cra i crC" Tcr . In contrast to a xy , which is only quantized 
if spin is conserved (in this case, the calculation of a xy and 
G s are completely equivalent, see section 3), the spin Chern 
number C s as defined in [62] remains quantized throughout 
the topological phase of the KM model with A R 7^ 
and/or disorder [62]. Whereas Sheng et al [62] suggest 
that (38) permits to distinguish two different quantum spin 
Hall phases (with the sign of a xy depending on the sign of 
the spin-orbit coupling), Fu and Kane [93] as well as Fukui 
and Hatsugai [178] argue that this distinction is merely a 
boundary effect, and that the state is fully characterized 
by a Z 2 index v — 1. The flux insertion related to the 
twisted boundary conditions can cause a spurious closing 
of the band gap on small systems [178, 179]. Although (35) 
and (39) are equivalent for the spin-conserved case, it is 
not clear if the above definition of the spin Chern number 
can also be applied (and remains quantized) to fractional 
topological insulators without spin conservation [138]. 

Equation (38) provides a rather general way to 
calculate the Z 2 index in correlated topological phases with 
time-reversal symmetry. However, its numerical evaluation 
can be quite involved or even technically impossible in the 
context of numerical approaches based on quantum Monte 
Carlo or quantum cluster methods. Consequently, a lot of 
effort has been devoted to developing simplified ways of 
calculating topological invariants. 

The starting point for one of the most fruitful recent 
developments is to express the Chern number in terms of 
the single-particle Green function [90, 171, 176, 180-182], 
dkdk 
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with k = oj (and cj real), fc = (k x , k y ), and = d/dk^ 
etc. G = G(iu), k) is the frequency and momentum- 
dependent single-particle Green function (in general a 
matrix in spin and orbital space), and G -1 its matrix 
inverse [171]. The trace is over the matrix indices of G 
and G _1 , and a summation over fj,, u, and p is implied. 
Equation (39) can be applied to spin-conserving quantum 
spin Hall systems by calculating C in each spin sector, 
see for example [72,75]. In this context, stochastic 
integration methods have been employed to sample the 
strongly varying integrands [183]. A simplification of (39), 
valid within the local approximation of dynamical mean- 
field theory, has been used to calculate a xy for the BHZH 
model [72], see figure 7. The Z 2 invariant of general, time- 
reversal invariant insulators can also be expressed in terms 
of the Green function; it is related to the electromagnetic 
polarization, and involves a five-dimensional integral [182]. 

Equation (39) can in principle be used directly to 
calculate the Chern number. Moreover, the formulation 
in terms of the single-particle Green function provides 
valuable insights into the bulk-boundary correspondence 
and interaction effects by analyzing the properties of 
G itself [90, 171, 184]. Most importantly, (39) implies 
that the Chern number can only change if a singularity 
appears in either G or G -1 [171, 184]. Considering first 
a noninteracting system, the relation G(iw,k) = [iuj — 
i/(k)] _1 reveals that for w on the real axis poles can only 
occur at to = 0, where they are related to zero eigenvalues 
of H and hence to zero-energy states. The change of C in 
this case is hence related to either a closing of the bulk band 
gap, or to an interface with a topologically different system 
where the zero-energy states correspond to edge states. In 
this way, the Green function can be used to derive the bulk- 
boundary correspondence principle [184]. 

In interacting systems, poles in the self-energy and 
hence in G _1 provide another mechanism for a change of 
the Chern number [107, 171], see (39). However, as pointed 
out in [171, 184], such Green function zeros are in general 
not related to a gap closing or to the existence of edge states. 
Consequently, the bulk-boundary correspondence can break 
down in the presence of interactions [171, 184]. Although 
intuitive, this picture based on the single-particle Green 
function may not capture interaction-generated collective 
excitations [171], such as gapless spinon excitations in a 
topological Mott insulator [29]. 

Whereas (39) is exact for noninteracting systems, 
its physical meaning is not yet fully understood in the 
interacting case. As discussed in [20, 171], for interacting 
systems, (39) relies on an adiabatic connection between 
the interacting state and a noninteracting, topological band 
insulator. It is therefore in general not applicable to 
interaction-driven states such as topological Mott insulators 
and fractional topological insulators. Although (39) 
remains a quantized, topological quantity even in the 
presence of interactions [171], its relation to physical 
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observables such as a s xy for general interacting states is an 
interesting open question [20, 171, 184]. On the other hand, 
(39) can be extended to disordered systems by replacing 
the physical momenta with phases representing twisted 
boundary conditions [20]. The Green function has also 
been used to study topological superconductors [185, 186], 
and correlated topological states in one [187] and three 
dimensions [188]. 

The multiple integrals, and the required knowledge of 
the Green function and its derivative, make the numerical 
evaluation of (39) rather difficult. Starting from (39), 
further simplifications have been achieved by Wang and 
Zhang [189] and by Wang et al [190], and the single- 
particle Green function again plays a key role. Wang and 
Zhang [189] showed that given an adiabatic connection 
to a noninteracting state (see also [20]), the topological 
invariant can be obtained from the Green function at uj = 
only [189], thereby simplifying practical calculations 
significantly. As shown in [189], the Chern number 
can be written in the form (13), with the Berry flux 
defined over the space of eigenvalues of the zero-frequency 
Green function. A related, physically more intuitive 
concept is the so-called topological Hamiltonian defined 
as h(k) = -G-^O,*;) = h (k) + £(0,fc) [175], 
where h a (k) is the noninteracting part of a many -body 
Hamiltonian, and E(0, k) is the zero-frequency part of the 
self-energy. Remarkably, provided there is an adiabatic 
connection to a noninteracting system, the topological 
Hamiltonian captures the topological invariant of the full, 
interacting problem [189]; see also [20, 175]. It therefore 
provides a highly nontrivial extension of the fact that for 
a noninteracting system, where — G _1 (0, k) = H(k), 
the Green function inherits the symmetries (time-reversal, 
particle -hole, and chiral symmetry) of the Hamiltonian, and 
can hence be used for a topological classification [184]. 
An illuminating discussion of why the zero-frequency self- 
energy is sufficient to calculate the topological invariant can 
be found in [175]. For practical applications, the power of 
the concept of the topological Hamiltonian is related to the 
possibility of calculating the Chern number by any method 
that is applicable to noninteracting systems. For example, 
given inversion symmetry, the Green function has to be 
determined only at the four time-reversal invariant momenta 
[190]. This method has been applied to the KMH model 
[107], and is quite similar in spirit to the work of Fu and 
Kane [7], in which the eigenvalues of the parity operator 
at the time-reversal invariant points in the Brillouin zone 
determine the topological invariant. 

In addition to this rather general approach, several 
other simplifications have been proposed and used in 
combination with numerical methods. The calculation of 
topological invariants becomes even simpler if the self- 
energy is local, i.e., independent of k. This approximation 
underlies the dynamical mean-field theory, and is expected 
to yield quite accurate results in high dimensions [191]. 



(a) (b) 
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Figure 19. The four possible adiabatic processes (pf(t), 4>_i(t) for 
inserting a n flux into a quantum spin Hall insulator [193]. Considering a 
closed loop around the flux, Faraday's law implies the existence of charge 
fluxons with charge ±e and spin 0, and spin fluxons with charge and 
spin ±1/2. (Reprinted with permission from [193]. Copyright 2008 by 
the American Physical Society). 

Wang et al [192] showed that a local self-energy that is 
diagonal in orbital space permits to write the Chern number 
as a product of the frequency-domain winding number and a 
Chern number of a mean-field Hamiltonian. The former is 
related to fluctuation-driven topological phase transitions, 
whereas the latter captures band structure effects. This 
scheme has been successfully applied to the KMH model 
[107]. For local self-energies that are not diagonal in 
orbital space, a simplification can be achieved using a pole 
expansion of the self-energy, and leads to an effective, 
noninteracting Hamiltonian whose topological invariant can 
be shown to be the same as that of the original, correlated 
model [183]. Applications of this method can be found 
in [73, 107, 183]. Interestingly, these calculations predict 
interaction-driven transitions to topological phases that are 
not captured by mean-field theory and hence related to 
quantum fluctuations [73, 107]. A possible drawback for 
applications to two-dimensional systems is the complete 
neglect of spatial correlations. 

Topological insulators can also be identified by their 
response to topological defects such as it fluxes [127, 
193], twisted boundary conditions [194], dislocations [69], 
vacancies [195], or disclinations [196]. Such defects, near 
which the topological properties are perturbed compared to 
the rest of the system, typically give rise to new, low-energy 
excitations with well-defined quantum numbers. Similar 
to the metallic edge states arising from open boundary 
conditions, the existence of these states is protected by 
time-reversal symmetry. They also have a close conceptual 
relation to the familiar soliton excitations in polyacetylene 
first proposed by Su et al [197]. For numerical methods, 
ir fluxes are of particular interest as they can be easily 
implemented, and will be the focus here. 

The existence of low-energy excitations as a result of 
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a it flux becomes quite apparent from the flux insertion 
argument by Qi and Zhang [193]; see also [12]. Consider 
a disk-shaped, noninteracting quantum spin Hall insulator 
that conserves spin [193]. Through a hole at the centre, 
a magnetic flux of magnitude hc/2e (equal to it in units 
where h = c = e = 1) is inserted adiabatically. This 
insertion can be described by the time-dependent fluxes 
(f> t (t), ^(t), with 4> a (t = 0) = and 4> a (t = 1) = ±tt. In 
total, there are four different possibilities to switch on the 
fluxes that differ by the signs of the fluxes 4>At) and <fi^(t), 
as shown in figure 19. Importantly, because the resulting 
Berry phases e 17r = e~ 17T = — 1, all four choices lead to 
the same final state [193]. Considering a loop around the 
flux tube, Faraday's law states that the flux gives rise to 
a tangential electric field, which for a quantum spin Hall 
insulator in turn causes a flow of charge perpendicular to 
the loop. The four different processes give rise to a doublet 
of states with charge Q — ±e and spin S z = (the so- 
called charge fluxons [127] or chargeons [193], figure 19(a) 
and (b)), and a pair of states with charge Q = and 
spin S z — ±1/2 (spin fluxons [127] or spinons [193], 
figure 19(c) and (d)). Given time-reversal and particle-hole 
symmetry, these states can be shown to lie at the middle of 
the band gap [127, 193]. Moreover, they are exponentially 
localized near the tt flux, and the two spinon states form 
a Kramers doublet related by time reversal, just like the 
two states of a spin 1/2. The whole argument can also be 
generalized to systems with only a Z 2 spin symmetry [193]. 

Fluxon states as a response to a it flux provide an 
alternative way to define the Z 2 topological invariant, and 
were suggested to be a useful tool to study correlated 
topological states in combination with numerical methods 
[127, 193]. In particular, there is no need to study systems 
with open boundaries, or to perform high-dimensional 
integrals over Green functions. On periodic lattices, it 
fluxes can be inserted in pairs, and are implemented in terms 
of Berry phase factors for hopping processes that cross the 
string or branch cut that connects two fluxes [101, 193, 
194]. A first study in this direction was recently carried 
out based on exact quantum Monte Carlo simulations of the 
KMH model [101]. In this work, it was established that the 
spinon states can be identified from a characteristic Curie 
law (given by N^/k B T for N fluxes) in the temperature- 
dependent magnetic susceptibility. This Curie law is 
related to the low-energy spinon doublet (in the presence 
of repulsive electronic interactions, charge fluxon states 
become gapped [127, 193]). As illustrated in figure 20, the 
presence or absence of a Curie law at low temperatures 
provides a clear distinction between a correlated quantum 
spin Hall state and an antiferromagnet with spontaneously 
broken time-reversal symmetry [101]. The presence of 
localized spinon excitations can also be visualized by 
means of the lattice-site dependent spin excitation spectrum 
[101]. Based on the flux insertion argument of figure 19, 
a generalization to fractional topological insulators seems 
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Figure 20. Magnetic susceptibility of the KMH model with one pair of tt 
fluxes from quantum Monte Carlo simulations onLxL lattices, (a) In the 
quantum spin Hall phase, midgap spinon states give rise to a Curie law at 
low temperatures. In contrast, such a Curie law is absent in the results for 
the magnetic insulating state shown in (b). Data taken from [101]. 



possible. 

Finally, quantum entanglement has in recent years 
emerged as one of the most useful tools to detect topological 
order in numerical calculations, most notably in the 
framework of the density matrix renormalization group and 
related approaches [198]. The entanglement entropy [199— 
201] is calculated as the partial trace over the density matrix 
of the ground state and provides a measure for the quantum 
correlations between different parts of the quantum system 
under consideration. It directly reveals the presence of 
excitations with exotic quantum numbers and can be used to 
identify, for example, quantum spin liquids [202, 203] and 
the exotic QSH* phase of the SI model [76]. Even more 
information is contained in the entanglement spectrum 
[204], corresponding to the spectrum of the reduced density 
matrix. In the context of topological insulators, it has been 
shown that topological phases possess protected gapless 
modes in their entanglement spectrum, which remain 
gapless even if the physical edge states are gapped out by 
interaction effects or perturbations that break time-reversal 
invariance [205,206]; for more details see [35,81] and 
references therein. Whereas the two-dimensional setting of 
quantum spin Hall insulators is typically difficult to study 
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with the density matrix renormalization group method, 
important progress has been made using tensor networks, 
and a suitable ansatz for fractional topological insulators 
has been proposed in [152]. 

7. Conclusions and outlook 

Despite their rather short history, time-reversal invariant 
topological insulators have already evolved into a multi- 
faceted field of research that is driven by interests ranging 
from future technological applications of the quantum 
spin Hall effect to the intricate physics of interaction- 
driven topological order. The present review has aimed at 
providing an overview of the fascinating topic of electronic 
correlation effects, which on their own have been at the 
focus of attention in the condensed matter community 
for decades. Similar to the fractional quantum Hall 
effect, the combination of electron-electron interactions 
and topological aspects in strongly-correlated topological 
insulators provides an ideal setting to explore beautiful new 
physics. 

Among the two-dimensional systems discussed here, 
the currently best understood correlated topological insu- 
lators are those with strong, intrinsic spin-orbit coupling 
that gives rise to a topological band insulator in the ab- 
sence of electronic correlations. If spin is conserved, such 
states have close conceptual relations to the integer quan- 
tum Hall effect. The impact of electronic interactions can in 
many cases be explored by means of sophisticated numeri- 
cal methods. Because time-reversal symmetry is at the heart 
of topological insulators, magnetic correlations induced by 
strong electron-electron interactions play a crucial role. For 
several of the available theoretical models, it is known that 
an extended topological insulator phase exists up to rather 
strong electronic interactions, and this phase appears to be 
adiabatically connected to a topological band insulator. At 
even stronger interactions, time-reversal symmetry can be 
spontaneously broken by the onset of long-range magnetic 
order. Generically, such a transition destroys the topolog- 
ical state but under special conditions, antiferromagnetic 
topological insulators with protected edge states exist. 

Correlation effects also play an important role for 
the helical edge states. Time-reversal symmetry provides 
protection against single-particle backscattering processes, 
but two-particle backscattering can give rise to long-range 
magnetic order at the edge, and to a breakdown of the 
bulk-boundary correspondence. However, whereas the 
relevant Luttinger parameter has been shown to become 
small enough for such a transition to occur, the necessary 
umklapp scattering is not allowed if spin is conserved. In 
the absence of magnetic impurities, the edges therefore 
retain their metallic character even for strong interactions. 
Helical edge states are also remarkably different from 
ordinary one-dimensional liquids concerning their stability 
with respect to magnetic impurities or disorder. 



Topological states that are not adiabatically connected 
to band insulators have been explored less thoroughly. 
Examples include topological Mott insulators, in which the 
topological state is generated from electronic correlations 
rather than from spin-orbit coupling, and fractional 
topological insulators that can, in the simplest case, be 
constructed by combining two fractional quantum Hall 
states in a time-reversal invariant way. In these cases, 
the resulting many-body state cannot be understood in the 
noninteracting limit, and intriguing phenomena such as 
quasiparticles with fractional charge and statistics appear. 
Moreover, while Rashba spin-orbit coupling — violating 
the conservation of total spin — is often considered as an 
experimental nuisance, theoretical results suggest that it 
may actually lead to exciting and not yet fully explored 
new physics. Topological Mott insulators and fractional 
topological insulators, as well as the general class of 
topologically ordered states with time-reversal symmetry, 
have so far been studied mostly using mean-field or 
topological field theory, and pose challenging problems 
such as finding suitable ways to determine whether or 
not a given state is topological. The identification 
and comprehensive numerical investigation of microscopic 
models is one of the future goals. Finally, the current search 
for experimental realizations, for example in cold-atom 
systems, will hopefully soon lead to a fruitful interplay of 
theory and experiment. 

Acknowledgments 

The authors are grateful to M. Bercx, B. Beri, J. Budich, 
N. Cooper, M. Daghofer, V. Gurarie, T. Lang, M. Laubach, 
Z. Meng, J. Moore, A. Muramatsu, S. Rachel, M. Schmidt, 
R. Thomale, B. Trauzettel, C. Varney, S. Wessel, and C. 
Xu for valuable discussions. Financial support from the 
DFG grants As 120/4-3 and Ho 4489/2-1 is gratefully 
acknowledged. 

References 

[1] K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 45, 494 
(1980). 

[2] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, 

Rev. Mod. Phys. 80, 1083 (2008). 
[3] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, 

Phys. Rev. Lett. 49, 405 (1982). 
[4] J. E. Avron, R. Seiler, and B. Simon, Phys. Rev. Lett. 51, 51 (1983). 
[5] M. Kohmoto, Annals of Physics 160, 343 (1985). 
[6] B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 
[7] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
[8] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
[9] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005). 
[10] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005). 
[11] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 106803 

(2007). 

[12] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (201 1). 
[13] H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Kleinman, and 
A. H. MacDonald, Phys. Rev. B 74, 165310 (2006). 



CONTENTS 



30 



[14] M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, L. W. 
Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 766 
(2007). 

[15] A. Roth, C. Briine, H. Buhmann, L. W. Molenkamp, J. Maciejko, 
X.-L. Qi, and S.-C. Zhang, Science 325, 294 (2009). 

[16] C. Briine, A. Roth, H. Buhmann, E. M. Hankiewicz, L. W. 
Molenkamp, J. Maciejko, X.-L. Qi, and S.-C. Zhang, Nat. Phys. 
8, 486 (2012). 

[17] B. A. Bernevig, T. L. Hughes, and S. Zhang, Science 314, 1757 
(2006). 

[18] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Phys. 

Rev. B 78, 195125 (2008). 
[19] A. Kitaev, AIP Conf. Proc. 1134, 22 (2009). 
[20] J. C. Budich and B. Trauzettel, arXiv: 1210.6672 (2012). 
[21] R.-J. Slager, A. Mesaros, V. Juricic, and J. Zaanen, 

arXiv: 1209.2610 (2012). 
[22] A. Shitade, H. Katsura, J. Kunes, X.-L. Qi, S.-C. Zhang, and N. 

Nagaosa, Phys. Rev. Lett. 102, 256403 (2009). 
[23] J. Reuther, R. Thomale, and S. Trebst, Phys. Rev. B 84, 100406 

(2011) . 

[24] J. Chaloupka, G. Jackeli, and G. Khaliullin, Phys. Rev. Lett. 105, 
027204 (2010). 

[25] J. Reuther, R. Thomale, and S. Rachel, Phys. Rev. B 86, 155127 

(2012) . 

[26] M. Kargarian, A. Langari, and G. A. Fiete, arXiv: 1207.2156 (2012). 
[27] L. Balents, Nature 464, 199 (2010). 

[28] S. Raghu, X. Qi, C. Honerkamp, and S. Zhang, Phys. Rev. Lett. 

100, 156401 (2008). 
[29] D. Pesin and L. Balents, Nat. Phys. 6, 376 (2010). 
[30] D. N. Sheng, Z.-C. Gu, K. Sun, and L. Sheng, Nat. Commun. 2, 

389 (2011). 

[31] J. W. F. Venderbos, S. Kourtis, J. van den Brink, and M. Daghofer, 

Phys. Rev. Lett. 108, 126405 (2012). 
[32] M. Levin and A. Stern, Phys. Rev. Lett. 103, 196803 (2009). 
[33] Y.-M. Lu and A. Vishwanath, Phys. Rev. B 86, 1251 19 (2012). 
[34] X.-G. Wen, Quantum Field Theory of Many-body Systems: From 

the Origin of Sound to an Origin of Light and Electrons (Oxford 

University Press, Oxford, 2004). 
[35] G. A. Fiete, V. Chua, X. Hu, M. Kargarian, R. Lundgren, A. Riiegg, 

J. Wen, and V. Zyuzin, Physica E 44, 845 (2012). 
[36] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 
[37] M. Z. Hasan and J. E. Moore, Annu. Rev. Cond. Mat. Phys. 2, 55 

(2011). 

[38] D. Xiao, W. Zhu, Y. Ran, N. Nagaosa, and S. Okamoto, Nat. 

Commun. 2,596 (2011). 
[39] A. Riiegg and G. A. Fiete, Phys. Rev. B 84, 201 103 (201 1). 
[40] R. Nandkishore and L. Levitov, Phys. Rev. B 82, 115124 (2010). 
[41] F. Zhang, J. Jung, G. A. Fiete, Q. Niu, and A. H. MacDonald, Phys. 

Rev. Lett. 106, 156801 (2011). 
[42] C. Weeks, J. Hu, J. Alicea, M. Franz, and R. Wu, Phys. Rev. X 1, 

021001 (2011). 

[43] K. K. Gomes, W. Mar, W. Ko, F. Guinea, and H. C. Manoharan, 

Nature 483, (2012). 
[44] J. Dalibard, F. Gerbier, G. Juzeliflnas, and P. Ohberg, Rev. Mod. 

Phys. 83, 1523 (2011). 
[45] B. Beri and N. R. Cooper, Phys. Rev. Lett. 107, 145301 (2011). 
[46] N. Goldman, I. Satija, P. Nikolic, A. Bermudez, M. A. Martin- 

Delgado, M. Lewenstein, and I. B. Spielman, Phys. Rev. Lett. 

105, 255302 (2010). 
[47] Y.-F. Wang, Z.-C. Gu, C.-D. Gong, and D. N. Sheng, Phys. Rev. 

Lett. 107, 146803 (2011). 
[48] J. Koch, A. A. Houck, K. L. Hur, and S. M. Girvin, Phys. Rev. A 

82,043811 (2010). 
[49] A. Petrescu, A. A. Houck, and K. L. Hur, arXiv:1206.1539 (2012). 
[50] J. E. Moore, Nature 464, 194 (2010). 

[51] M. Dzero, K. Sun, V. Galitski, and P. Coleman, Phys. Rev. Lett. 

104, 106408 (2010). 
[52] X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Phys. Rev. B 74, 085308 

(2006). 

[53] X. Qi, T. L. Hughes, and S. Zhang, Phys. Rev. B 78, 195424 (2008). 



[54] C. Wu, B. A. Bernevig, and S.-C. Zhang, Phys. Rev. Lett. 96, 
106401 (2006). 

[55] C. Xu and J. E. Moore, Phys. Rev. B 73, 045322 (2006). 
[56] A. M. Essin and V. Gurarie, Phys. Rev. B 85, 195116 (2012). 
[57] Z. Y. Meng, T. C. Lang, S. Wessel, F. F. Assaad, and A. Muramatsu, 

Nature 464, 847 (2010). 
[58] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, 

S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 

666 (2004). 

[59] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and 

A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 
[60] J. Hubbard, Proc. R. Soc. London 276, 238 (1963). 
[61] P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[62] D. N. Sheng, Z. Y. Weng, L. Sheng, and F. D. M. Haldane, Phys. 

Rev. Lett. 97, 036808 (2006). 
[63] S. Rachel and K. Le Hur, Phys. Rev. B 82, 075106 (2010). 
[64] K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Phys. Rev. Lett. 

103, 046811 (2009). 
[65] F. D. M. Haldane, Phys. Rev. Lett. 61, 1029 (1988). 
[66] L. Wang, H. Shi, S. Zhang, X. Wang, X. Dai, and X. C. Xie, 

arXiv:1012.5163 (2010). 
[67] C. N. Varney, K. Sun, M. Rigol, and V. Galitski, Phys. Rev. B 82, 

115125 (2010). 

[68] C. N. Varney, K. Sun, M. Rigol, and V. Galitski, Phys. Rev. B 84, 
241105 (2011). 

[69] V. Juricic, A. Mesaros, R.-J. Slager, and J. Zaanen, Phys. Rev. Lett. 

108, 106403 (2012). 
[70] A. Medhi, V. B. Shenoy, and H. R. Krishnamurthy, Phys. Rev. B 85, 

235449 (2012). 

[71] T. L. Schmidt, S. Rachel, F. von Oppen, and L. I. Glazman, Phys. 

Rev. Lett. 108, 156402 (2012). 
[72] T. Yoshida, S. Fujimoto, and N. Kawakami, Phys. Rev. B 85, 

125113 (2012). 

[73] L. Wang, X. Dai, and X. C. Xie, Eur. Phys. Lett. 98, 57001 (2012). 
[74] Y. Tada, R. Peters, M. Oshikawa, A. Koga, N. Kawakami, and S. 

Fujimoto, Phys. Rev. B 85, 165138 (2012). 
[75] T. Yoshida, R. Peters, S. Fujimoto, and N. Kawakami, 

arXiv: 1207 .4547 (2012). 
[76] A. Riiegg and G. A. Fiete, Phys. Rev. Lett. 108, 046401 (2012). 
[77] H.-M. Guo and M. Franz, Phys. Rev. B 80, 1 13 102 (2009). 
[78] A. Riiegg, J. Wen, and G. A. Fiete, Phys. Rev. B 81, 205 1 15 (2010). 
[79] C. Weeks and M. Franz, Phys. Rev. B 82, 085310 (2010). 
[80] M. Kargarian and G. A. Fiete, Phys. Rev. B 82, 085106 (2010). 
[81] X. Hu, M. Kargarian, and G. A. Fiete, Phys. Rev. B 84, 155116 

(2011). 

[82] A. Kitaev, Annals of Physics 303, 2 (2003). 
[83] A. Kitaev, Annals of Physics 321, 2 (2006). 

[84] M. A. N. Araujo, E. V. Castro, and P. D. Sacramento, 

arXiv:1208.1289 (2012). 
[85] J. Goryo and N. Maeda, J. Phys. Soc. Jpn. 80, 044707 (201 1). 
[86] D. Cocks, P. P. Orth, S. Rachel, M. Buchhold, K. Le Hur, and W. 

Hofstetter, arXiv:1204.4171 (2012). 
[87] Q. Niu, D. J. Thouless, and Y.-S. Wu, Phys. Rev. B 31, 3372 (1985). 
[88] D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 1959 

(2010). 

[89] W.-Y. Hsiang and D.-H. Lee, Phys. Rev. A 64, 052101 (2001). 
[90] G. E. Volovik, The Universe in a Helium Droplet (Clarendon Press, 
Oxford, 2003). 

[91] D. N. Sheng, L. Balents, and Z. Wang, Phys. Rev. Lett. 91, 1 16802 
(2003). 

[92] G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 
[93] L. Fu and C. L. Kane, Phys. Rev. B 74, 195312 (2006). 
[94] X. G. Wen, Phys. Rev. B 41, 12838 (1990). 
[95] X. G. Wen, Int. J. Mod. Phys. B 6, 1711 (1992). 
[96] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 (1992). 
[97] M. Onoda, Y. Avishai, and N. Nagaosa, Phys. Rev. Lett. 98, 076802 
(2007). 

[98] T. Neupert, L. Santos, S. Ryu, C. Chamon, and C. Mudry, Phys. 

Rev. B 84, 165107 (2011). 
[99] M. Hohenadler, T. C. Lang, and F. F. Assaad, Phys. Rev. Lett. 106, 



CONTENTS 



31 



100403 (201 1), erratum (in press). 
100] M. Hohenadler, Z. Y. Meng, T. C. Lang, S. Wessel, A. Muramatsu, 

and F. F. Assaad, Phys. Rev. B 85, 115132 (2012). 
101] F. F. Assaad, M. Bercx, and M. Hohenadler, arXiv: 1204.4728 

(2012). 

102] C. Griset and C. Xu, Phys. Rev. B 85, 045123 (2012). 
103] D. Zheng, G.-M. Zhang, and C. Wu, Phys. Rev. B 84, 205121 
(2011). 

104] Y. Yamaji and M. Imada, Phys. Rev. B 83, 205122 (201 1). 
105] W. Wu, S. Rachel, W.-M. Liu, and K. Le Hur, Phys. Rev. B 85, 
205102 (2012). 

106] S.-L. Yu, X. C. Xie, and J.-X. Li, Phys. Rev. Lett. 107, 010401 

(2011) . 

107] J. C. Budich, R. Thomale, G. Li, M. Laubach, and S.-C. Zhang, 

arXiv: 1203.2928 (2012). 
108] J. Wen, M. Kargarian, A. Vaezi, and G. A. Fiete, Phys. Rev. B 84, 

235149 (2011). 

109] M. Mardani, M.-S. Vaezi, and A. Vaezi, arXiv: 1111 .5980 (201 1). 
1 10] A. Vaezi, M. Mashkoori, and M. Hosseini, Phys. Rev. B 85, 195126 

(2012) . 

1 1 1] D.-H. Lee, Phys. Rev. Lett. 107, 166806 (201 1). 

1 12] S. Sorella, Y. Otsuka, and S. Yunoki, arXiv: 1207. 1783 (2012). 

113] M. Hohenadler and F. F. Assaad, Phys. Rev. B 85, 081 106 (2012), 

erratum (in press). 
114] R. S. K. Mong, A. M. Essin, and J. E. Moore, Phys. Rev. B 81, 

245209 (2010). 

115] M. Campostrini, M. Hasenbusch, A. Pelissetto, and E. Vicari, Phys. 

Rev. B 74, 144506 (2006). 
116] T. Grover and A. Vishwanath, arXiv:1206.1332 (2012). 
117] V. Anisimov, I. Nekrasov, D. Kondakov, T. Rice, and M. Sigrist, 

Eur. Phys. J. B 25, 191 (2002). 
118] A. Liebsch, Phys. Rev. Lett. 91, 226401 (2003). 
119] A. Koga, N. Kawakami, T. M. Rice, and M. Sigrist, Phys. Rev. Lett. 

92, 216402 (2004). 
120] R. Arita and K. Held, Phys. Rev. B 72, 201 102 (2005). 
121] H. Guo, S. Feng, and S.-Q. Shen, Phys. Rev. B 83, 045114 (2011). 
122] J. He, Y.-H. Zong, S.-P. Kou, Y. Liang, and S. Feng, Phys. Rev. B 

84,035127(2011). 
123] J. He, B. Wang, and S.-P. Kou, arXiv:1204.4766 (2012). 
124] T. Neupert, L. Santos, S. Ryu, C. Chamon, and C. Mudry, Phys. 

Rev. Lett. 108, 046806 (2012). 
125] I. F. Herbut, Phys. Rev. Lett. 97, 146401 (2006). 
126] C. Honerkamp, Phys. Rev. Lett. 100, 146404 (2008). 
127] Y. Ran, A. Vishwanath, and D.-H. Lee, Phys. Rev. Lett. 101, 

086801 (2008). 

128] R. Nandkishore, M. A. Metlitski, and T. Senthil, Phys. Rev. B 86, 
045128 (2012). 

129] C. Wu and S.-C. Zhang, Phys. Rev. Lett. 93, 036403 (2004). 
130] C. Weeks and M. Franz, Phys. Rev. B 81, 085105 (2010). 
131] Q. Liu, H. Yao, and T. Ma, Phys. Rev. B 82, 045102 (2010). 
132] J. Wen, A. Riiegg, C.-C. J. Wang, and G. A. Fiete, Phys. Rev. B 82, 
075125 (2010). 

133] Y. Zhang, Y. Ran, and A. Vishwanath, Phys. Rev. B 79, 245331 
(2009). 

134] M. W. Young, S.-S. Lee, and C. Kallin, Phys. Rev. B 78, 125316 
(2008). 

135] B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 106802 
(2006). 

136] E. Tang, J.-W. Mei, and X.-G. Wen, Phys. Rev. Lett. 106, 236802 
(2011). 

137] K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, Phys. Rev. Lett. 106, 
236803 (2011). 

138] T. Neupert, L. Santos, C. Chamon, and C. Mudry, Phys. Rev. Lett. 

106,236804(2011). 
139] M. Levin, F. J. Burnell, M. Koch-Janusz, and A. Stern, Phys. Rev. 

B 84,235145 (2011). 
140] A. M. Lauchli, Z. Liu, E. J. Berholtz, and R. Moessner, 

arXiv: 1207.6094 (2012). 
141] Y.-H. Wu, J. K. Jain, and K. Sun, Phys. Rev. B 86, 165129 (2012). 
142] G. Y. Cho and J. E. Moore, Annals of Physics 326, 1515 (201 1). 



[143] M. Levin and A. Stern, Phys. Rev. B 86, 115131 (2012). 
[144] Y.-M. Lu and Y. Ran, Phys. Rev. B 85, 165134 (2012). 
[145] D. Ferraro and G. Viola, arXiv:l 1 12.5399 . 
[146] R. Roy, arXiv: 1208.2055 (2012). 
[147] P. Nikolic, arXiv: 1108.5388 (2011). 
[148] P. Nikolic, arXiv:1206.1055 (2012). 

[149] A. G. Grushin, T. Neupert, C. Chamon, and C. Mudry, 

arXiv: 1207.4097 (2012). 
[150] B. Beri and N. R. Cooper, Phys. Rev. Lett. 108, 206804 (2012). 
[151] L. Tianhan, C. Repellin, B. A. Bernevig, and N. Regnault, 

arXiv: 1206.2626 (2012). 
[152] B. Beri and N. R. Cooper, Phys. Rev. Lett. 106, 156401 (201 1). 
[153] J. Voit, Rep. Prog. Phy. 57, 977 (94). 

[154] G. Tkachov and E. M. Hankiewicz, arXiv: 1208. 1466 (2012). 
[155] G. Autes and O. V. Yazyev, arXiv:1210.4818 (2012). 
[156] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Nat. Phys. 4, 273 (2008). 
[157] J. C. Budich, F. Dolcini, P. Recher, and B. Trauzettel, Phys. Rev. 

Lett. 108, 086602 (2012). 
[158] A. Medhi and V. B. Shenoy, J. Phys: Condens. Matter 24, 355001 

(2012). 

[159] A. V. Moroz, K. V. Samokhin, and C. H. W. Barnes, Phys. Rev. Lett. 

84, 4164 (2000). 
[160] A. Iucci, Phys. Rev. B 68, 075107 (2003). 

[161] Y. A. Bychkov and E. I. Rashba, J. Phys. C: Solid State Phys. 17, 
6039 (1984). 

[162] J. Maciejko, C. Liu, Y. Oreg, X.-L. Qi, C. Wu, and S.-C. Zhang, 

Phys. Rev. Lett. 102, 256803 (2009). 
[163] J. Maciejko, Phys. Rev. B 85, 245108 (2012). 
[164] F. Crepin, J. C. Budich, F. Dolcini, P. Recher, and B. Trauzettel, 

Phys. Rev. B 86, 121 106 (2012). 
[165] D. Soriano and J. Fernandez-Rossier, Phys. Rev. B 82, 161302 

(2010). 

[166] H. Feldner, Z. Y. Meng, T. C. Lang, F. F. Assaad, S. Wessel, and A. 

Honecker, Phys. Rev. Lett. 106, 226401 (2011). 
[167] D. J. Luitz, F. F. Assaad, and M. J. Schmidt, Phys. Rev. B 83, 

195432 (2011). 
[168] M. J. Schmidt, arXiv:1207.3801 (2012). 

[169] B. Braunecker, C. Bena, and P. Simon, Phys. Rev. B 85, 035136 
(2012). 

[170] X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Phys. Rev. B 74, 045125 
(2006). 

[171] V. Gurarie, Phys. Rev. B 83, 085426 (201 1). 

[172] A. Furusaki and N. Nagaosa, Phys. Rev. Lett. 72, 892 (1994). 

[173] J. Yuan, J.-H. Gao, W.-Q. Chen, F. Ye, Y. Zhou, and F.-C. Zhang, 

Phys. Rev. B 86, 104505 (2012). 
[174] R. R. P. Singh and R. T. Scalettar, Phys. Rev. Lett. 66, 3203 (1991). 
[175] Z. Wang and B. Yan, arXiv: 1207.7341 (2012). 
[176] K. Ishikawa and T. Matsuyama, Nucl. Phys. B 280, 523 (1987). 
[177] S. Kourtis, J. W. F. Venderbos, and M. Daghofer, arXiv: 1208.3481 

(2012). 

[178] T. Fukui and Y. Hatsugai, Phys. Rev. B 75, 121403 (2007). 
[179] S.-S. Lee and S. Ryu, Phys. Rev. Lett. 100, 186807 (2008). 
[180] K. Ishikawa and T. Matsuyama, Z. Physik C 33, 41 (1986). 
[181] H. So, Prog. Theor. Phys. 74, 585 (1985). 

[182] Z. Wang, X.-L. Qi, and S.-C. Zhang, Phys. Rev. Lett. 105, 256803 
(2010). 

[183] L. Wang, H. Jiang, X. Dai, and X. C. Xie, Phys. Rev. B 85, 235135 
(2012). 

[184] A. M. Essin and V. Gurarie, Phys. Rev. B 84, 125132 (2011). 
[185] Z. Wang and S.-C. Zhang, Phys. Rev. B 86, 165116 (2012). 
[186] J. C. Budich and B. Trauzettel, arXiv:1207.1104 (2012). 
[187] S. R. Manmana, A. M. Essin, R. M. Noack, and V. Gurarie, 

arXiv: 1205.5095 (2012). 
[188] A. Go, W. Witczak-Krempa, G. S. Jeon, K. Park, and Y. B. Kim, 

Phys. Rev. Lett. 109, 066401 (2012). 
[189] Z. Wang and S.-C. Zhang, Phys. Rev. X 2, 031008 (2012). 
[190] Z. Wang, X.-L. Qi, and S.-C. Zhang, Phys. Rev. B 85, 165126 

(2012). 

[191] A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. Mod. 
Phys. 68, 13 (1996). 



CONTENTS 



[192] L. Wang, X. Dai, and X. C. Xie, Phys. Rev. B 84, 205116 (2011). 
[193] X.-L. Qi and S.-C. Zhang, Phys. Rev. Lett. 101, 086802 (2008). 
[194] D.-H. Lee, G.-M. Zhang, and T. Xiang, Phys. Rev. Lett. 99, 196805 
(2007). 

[195] J. He, Y.-X. Zhu, Y.-J. Wu, L.-F. Liu, Y. Liang, and S.-P. Kou, 

arXiv: 1210.0266 (2012). 
[196] A. Riiegg and C. Lin, arXiv: 1208.3023 (2012). 
[197] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 42, 

1698 (1979). 

[198] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 
[199] G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90, 
227902 (2003). 

[200] A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 1 10404 (2006). 
[201] M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 (2006). 
[202] S. Depenbrock, I. P. McCulloch, and U. Schollwock, Phys. Rev. 

Lett. 109, 067201 (2012). 
[203] H.-C. Jiang, Z. Wang, and L. Balents, arXiv:1205.4289 (2012). 
[204] H. Li and F. D. M. Haldane, Phys. Rev. Lett. 101, 010504 (2008). 
[205] A. M. Turner, Y. Zhang, and A. Vishwanath, Phys. Rev. B 82, 

241102 (2010). 

[206] X.-L. Qi, H. Katsura, and A. W. W. Ludwig, Phys. Rev. Lett. 108, 
196402 (2012). 



